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INTRODUCTION 


Researchers at the Langley Research Center are currently developing and applying 
the latest optimization capability to the multidisciplinary optimization of aircraft 
and spacecraft (ref. 1). A key part of this effort is to efficiently calculate 
structural-sensitivity derivatives which quantify the change in a behavior variable 
with respect to a structural parameter and which are used in the following applica- 
tions: (1) to act as input to optimization algorithms; (2) to enhance response anal- 

ysis programs which aid engineering judgment leading to design modifications? (3) to 
guide the modification of a finite-element model to better correlate analytical and 
test results; and (4) to approximate structural response by using Taylor series 
expansions. 

The most basic and straightforward approach to sensitivity analysis is the 
finite-difference method; however, it is computationally slow and its accuracy must 
be verified by convergence checks. Analytical methods for calculating exact deriva- 
tives from the governing equations have been developed (refs. 2 to 10). These 
methods greatly reduce the computational effort but are somewhat cumbersome to imple- 
ment for bending-type elements (ref. 4). Most recently, a semianalytical method for 
calculating derivatives has been developed which has the generality and programing 
ease of the finite-difference method while retaining much of the efficiency of the 
analytical method (ref. 11). 

The purpose of this paper is to describe the implementation and verification of 
the aforementioned methods for computing sensitivity derivatives in the structural 
finite-element computer program currently used in multidisciplinary optimization 
studies. The program, denoted the Engineering Analysis Language (EAL) System 
(ref. 12), is similar to its predecessor, SPAR (ref. 13). It is a modular system of 
individual analysis processors which may be used in any appropriate sequence to 
perform a variety of analyses. The EAL System differs from SPAR by providing 
FORTRAN-like commands which permit branching, testing data, looping, and calling 
runstreams (similar to calling FORTRAN subroutines). These capabilities permit the 
implementation of sensitivity calculations without changing the basic program or 
requiring user-written subroutines in separate programs (as in ref. 4) • Further, use 
of. the EAL System avoids the need for extensive operating system control commands as 
used previously (ref. 4) and thus assures machine independence of the resulting 
system. 

This paper draws on results from references 4, 8, and 11 for the basic method- 
ology and presents EAL input runstreams which calculate derivatives of displacements, 
stresses, vibration frequencies and mode shapes, and buckling loads and mode shapes 
with respect to structural variables in the model. The variables are sectional 
properties including thicknesses, cross-sectional areas, and moments of inertia. 
Results are presented and comparisons are made among analytical, semianalytical, and 
finite-difference methods for the following four structural conf igurations: a swept 

wing, a box beam, a stiffened cylinder with a cutout, and a space radiometer-antenna 
truss. 



SYMBOLS AND ABBREVIATIONS 


A cross-sectional area 

E modulus of elasticity 

{f} pseudo load vector 

[F] matrix of pseudo load vectors (eq. (21)) 

f applied load 

{f} applied load vector 

[G] matrix which relates stresses and temperatures (eq. (26) ) 

[I] identity matrix 

1^,1 2 moments of inertia 

J Q polar moment of inertia 

[K] stiffness matrix 

[Kg] geometric stiffness matrix 

[M] mass matrix 

Pj structural definition parameter 

{Q} particular solution of equation (35) 

r independent design variable for tube element (eq. (40) ) 

r^,r Q inner and outer radii of tube element 

[S] stress-displacement matrix 

T temperature 

T q stress-free temperature 

t^ element thickness 

u,v,w displacements in x, y, and z directions, respectively 
{u} displacement vector 

{v} vector of independent design variables 

{v} vector of dependent design variables 

x,y,z Cartesian coordinates 

vector of linear thermal expansion coefficients 
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vectors of buckling loads and vibration frequencies, respectively 
v Poisson’s ratio 

[$]/[<!>] matrices of buckling and vibration mode shapes, respectively 

p density 

{cr} stress vector 

Superscripts: 

n new (perturbed) 

o original (unperturbed) 

Abbreviations: 

NDDV number of dependent design variables 

NIDV number of independent design variables (NDV in computer listing) 

NPOL degree of polynominal used in linking (eq. (1 ) ) 

NSDP number of structural definition parameters 

DESCRIPTION OF THE ENGINEERING ANALYSIS LANGUAGE (EAL) SYSTEM 

The EAL finite-element analysis system (ref. 12) evolved from the SPAR computer 
program (ref. 13). As indicated in figure 1(a), the EAL System contains individual 
processors which communicate through a data base consisting of one or more libraries 
of data sets. The data sets typically contain data describing the finite-element 
model of the structure (e.g. , node point coordinates and material properties) as well 
as response information such as displacements and stresses. All data base communi- 
cations between processors are in terms of data sets. A set of data-handling util- 
ities transfers data between the processors in central memory of the computer and the 
data base on auxiliary storage (ref. 14). Contents of individual data sets for SPAR 
are compiled and listed in reference 15. (In some instances these data sets may be 
slightly different from corresponding EAL data sets.) 

A list of the EAL processors is shown in the uppermost block of figure 1(a). 

The functions of the various processors are described in table 1. The processors may 
be executed in any appropriate sequence. A sequence of processor executions is 
denoted a runstream and may be defined and stored in the data base as a runstream 
data set. Runstreams and runstream data sets may be nested within an input file. 

The EAL system differs from SPAR in its use of a set of flexible FORTRAN-like 
statements, denoted executive control system (ECS) commands, which allow branching, 
testing data, looping, and calling runstreams (similar to calling FORTRAN sub- 
routines) . The ECS commands are always preceded by an asterisk (*) and are used to 
execute processors; for example, the ECS command to execute the TAB processor is 
*XQT TAB. A list of the ECS commands and their functions is given in table 8-1 of 
reference 12. The ECS commands are also used to test variables (denoted registers). 
Registers may be input data, output data, or simple variables defined and manipulated 
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by the EAL user. If the value of a register is regarded as a variable which changes 
during an EAL analysis , the current value of the register is defined as a surrogate. 

A surrogate is always enclosed in quotation marks. Registers and surrogates are 
defined and manipulated by the use of register action commands (RAC). A RAC is 
always preceded by an exclamation point (!) in the EAL input. For example, 

! A=ABS ( "STR" ) assigns the absolute value of the surrogate "STR" to register A. A 
complete list of the register action commands is given in table 7.3-1 of 
reference 1 2. 

The makeup of a sample EAL input file is illustrated in figure 1(b). The file 
begins with an ECS command to execute processor TAB using appropriate input data 
following the command? second, there is a call to execute runstream data set CHNG DV? 
third, there is execution of processor ELD using appropriate input data; and finally, 
there is execution of processors of K and INV, neither of which require user input. 

The makeup of the sample runstream data set CHNG DV is shown in figure 1(c). 

The function of this runstream is to loop over the membrane elements in the model 
(looping is controlled by the JLZ command), extract the membrane thicknesses from 
data set DV DFN and store them in register DV (using the DS register action command), 
and finally update that part of the input file in which the membrane thicknesses are 
defined. The updated portion of the file follows the SA identifier and contains the 
membrane element numbers ("NI" ) and the corresponding thicknesses ("DV" ). Following 
the satisfaction of the loop and the arrival at LABEL 2, control returns to the call- 
ing runstream by means of the RETURN command. 


SCOPE OF CAPABILITY FOR CALCULATING STRUCTURAL-SENSITIVITY DERIVATIVES 
Methods Used and Overall Capability 

In the present work, structural behavior quantities include displacements {u}, 
stresses {a}, buckling loads and mode shapes {X} and [3>], and vibration fre- 
quencies and mode shapes {co 2 } and [<J> ] . Structural variables include element 
thicknesses, cross-sectional areas, and moments of inertia. Ihe following methods 
for calculating derivatives of structural behavior quantities have been implemented: 
(1) an analytical method, (2) a semi analytical method (also called the indirect 
method in refs. 11 and 16), and (3) a finite -difference method (also called the 
direct method in ref. 11). A list of applicable EAL finite elements is given in 
table 2. The analytical method was not implemented for the elements which include 
bending deformation (E32, E33, E42, and E43) because of the algebraic complexity 
(discussed in the following section). 


Design-Variable Definition 

Structural modifications are specified as changes to certain structural quanti- 
ties called design variables, which are related to section properties or mass proper- 
ties of the finite-element model. A structural definition parameter (ref. 4) is 
defined to be a parameter which has a linear contribution to the stiffness matrix or 
the mass matrix, or both, of individual finite elements in the structural model. The 
design variables can be identical to or have a one-to-one relationship to the struc- 
tural definition parameters. For example, the cross-sectional area of a rod element 
is a structural definition parameter, and the areas of several elements in the struc- 
ture could be equal to a single design variable. In some instances, there is a non- 
linear relationship between the structural definition parameters and the design 
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variables, for example, when the moment of inertia per unit width of a plate 
I 1 = t /I 2 is a structural definition parameter and the place thickness t is a 
design variable. Finally, in some optimization techniques (e.g., ref. 3), design- 
variable linking is used to reduce the number of independent design variables. In 
this instance independent design variables can be linearly or nonlinearly related to 
the dependent design variables, which in turn can be linearly or nonlinearly related 
to the structural definition parameters. Chain-rule partial differentiation is then 
used to compute the required derivatives (ref. 4). 


Design-Variable Linking 

When linking is used, the number of dependent design variables (NDDV) is larger 
than the number of independent design variables (NIDV) , and the two are related via a 
mathematical relationship which is often linear (ref. 17). In the present work the 
following nonlinear design-variable linking algorithm is used: 


v. = C . 
1 oi 


NPOL NIDV 

+ L E c 


, V, 


m 


m=l k=1 "k'i k 


(i =1, 2, ..., NDDV) 


( 1 ) 


where NPOL is the degree of the polynomial expression for each dependent design vari- 
able, is the ith linking coefficient, C oi is the ith additive constant, 

v^ is the ith dependent design variable, and V k is the independent design 
variable. With matrix notation, equation (1) may be written 


NPOL 

{v} = {c} + y [C] {v} m 

o ' m 

m=1 


( 2 ) 


STIFFNESS AND MASS MATRIX DERIVATIVES 

Ihe derivatives of the stiffness matrix [K] and of the mass matrix [M] with 
respect to the independent design variable V^ are needed in the subsequent calcu- 
lations of sensitivity derivatives. Two methods are used for calculating these 
matrix derivatives: analytical and finite differences. 


Analytical Method 

When the stiffness or the mass matrix of a finite element is linearly related to 
a design variable, the analytical derivative of the matrix is obtained by simply set- 
ting the design variable to unity and calculating the corresponding matrix. When a 
nonlinear relationship exists between the stiffness or the mass matrix and the design 
variables or when design -variable linking is used, or both, the situation is more 
complicated and chain-rule partial differentiation is needed to formulate generalized 
analytical expressions. For a beam element the stiffness matrix is linearly related 
to four structural definition parameters (A, 1 ^, I 2 , and J q ) . Analytical 
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expressions may be written to relate the design variables to these structural 
definition parameters, and the derivative of the stiffness matrix can be calculated 
by using 


d[K] _ 5[K] 5A d[K] 5I 1 5 [K] 5I 2 5[K] SJ Q 

dv. 5A dv. + dl 5v. + dl ‘5v7 + dJ dv. 

i i 1 i 2 i o i 


or 


S[K] 

8v, 


NSDP 


Z 


aw dP j 

BP. Sv, 
3 i 


(3) 


(4) 


where NSDP is the number of structural definition parameters, 3[K]/5Pj is computed 
by substituting unity for P j , and dPj/dv^ is calculated by differentiating the 
analytical expression relating v^ to Pj. The analogous expression for mass matrix 
derivatives is 


d[M] 

5v. 


NSDP 


E 


5p 

B[M] j 

dP. 

D 


Bv. 


(5) 


If design-variable linking is used, the independent and dependent design vari- 
ables are related by equation (1) and the expressions for derivatives with respect to 
the independent design variables are 


and 


BIK] 

5V 


8[M] 

dv. 


NDDV NSDP 



NDDV NSDP 



dW 

5P. 

3 

Bv. 

i 

BP. 

3 

Bv. 

i 

5 \ 


d[M] 

o> 

j* 

Bv. 

i 

BP. 

Bv. 

dv 

3 


k 


( 6 ) 


(7) 


Example of Analytical Derivative of Stiffness Matrix 

To crystallize the ideas and terms involved, differentiation of the stiffness 
matrix for a beam modeled by three E23 channel-section elements (fig. 2) is per- 
formed. The structural definition parameters consist of the cross-sectional area A 
and the moments of inertia 1 ^, I 2 , and J Q . The dependent design variables are the 

dimensions of the channel , B 2 , and t. There are three independent design 
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variables, each associated with an element and denoted , V 2 , and V^. The 
dependent design variables are linked to the independent design variables by a 
9 x 3 matrix [C] . Thus, for the structure in figure 2, 


{p} = {A I I J } 

1 z o 


(8) 


{v } 1 = ^ B 1 B 2 t 


B iV 


B 1 B 2 t} 


(9) 


Element 1 Element 2 Element 3 


{v } 1 = {v 1 v 2 v 3 > 


{v} = [C]{V} 


( 10 ) 

( 11 ) 


The structural definition parameters are related to the dependent design variables by 
equations given in reference 4, as follows: 




P- = 


A = (2B i + B 2 )t 


P„ = 




h = 


B 1 (B 2 + 2 t ) 3 (B 1 - t)B 2 3 


12 


12 
x 2 


3 2 3 

2tB /B \ B t J 

h ’ — + 2tB 1— - c + — * B 2 tl = 




(1 2 ) 


P. = J = 4-( 2B i + B o ) fc3 
4 o 3 1 2 


c = 


B i 2t 


j 


From equation ( 6 ), 


9 4 


5[K] 

dv. 



i =1 j =1 


Ar „, ap. 5v. 
OIK] i i 


5P. dv. 5V 
3 l k 


(13) 
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where 


8[K] 

6P. 

D 


[K] 


P.=1 

3 


(j = 1, 2, 3 , 4) 


(14) 


aPj/avi is obtained from differentiating equations (12) and noting equation (9), and 
5 v i /3V k is obtained from differentiating equation (11). 


Finite-Difference Method 

The stiffness and mass matrix derivatives can also be computed by finite 
differences, as follows: 


5 [K] _ A[K] _ 


av. 


Av 


[K] k n - [K]° 

n o 

V. - V, 
k k 


(15) 


a[M] _ a[m] 
av “ av 

k k 





(16) 


where [K] k n and [M] k n are P ertur ked stiffness and mass matrices (formed by 
incrementing the kth independent design variable), [K]° and [M]° are the original 
(unperturbed) stiffness and mass matrices, and V^ n and V k ° are the perturbed and 
original values of the kth independent design variable. Since this method operates 
with the independent design variables directly, it does not require looping over the 
dependent design variables or the structural parameters indicated by equations (6) 
and (7). 


Finite Differences With LSK Processor 

In some instances the finite-difference calculation of the stiffness matrix may 
be improved by use of the LSK processor (table 1). This processor selects appro- 
priate combinations of elements (submatrices) of the global stiffness matrix [K] in 
a sparse matrix form (called LS-format) and performs the necessary finite differenc- 
ing. For this case the differencing must be with respect to the dependent design 
variables, as follows: 


&[K] s A[K] 

5 v . Av . 

i i 


[K]. n 


n 

v. 

i 


- [K] 


o 


- v. 
1 


o 


(17) 
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If design-variable linking is used (eg. (1)), 


a[Ki 


dv 


K 


NDDV 




5v. 3v 
1 k 


( 18 ) 


where 3v i /9V k is obtained by differentiating equation ( 1 ) • Using the LSK processor 
to operate on selected portions of the global stiffness matrix is more efficient than 
summing and multiplying stiffness matrices having the dimensions of the entire struc- 
ture. Hence, depending on the ratio of the number of dependent design variables 
(NDDV) to the number of independent design variables (NIDV), the finite-difference 
methods of equation (15) or of equations (17) and (18) may be more efficient than the 
analytical method of equation (6). A large ratio of dependent to independent design 
variables and a large number of structural parameters needed to specify a particular 
element favor the use of one of the finite-difference methods for calculating 
3[K]/dV k . 


DISPLACEMENT DERIVATIVES 
Analytical Method 

In the analytical method, derivatives are computed from the governing finite- 
element equations. For static finite-element structural analysis, the equilibrium 
equation is 


[K] {u} = {f} 


(19) 


where {u} is the vector of nodal displacements and {f} is the applied load vec- 
tor. Differentiating equation (19) with respect to the independent design variable 
V k gives 


or 


[K] 


5{u} _ 5{f } d[K] 


dV, 


av. 


av. 


{u} = {F> 



[F] 


( 20 ) 


( 21 ) 


where {f}^ is the kth pseudo applied load vector and [5u/dV] and [F] are, 
respectively, matrices for which the columns are the displacement derivatives and 
pseudo load vectors which correspond to individual independent design variables. If 
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the applied load vector {f} is not a function of the design variables, then the 
term 9{f}/9V^ is equal to zero. ^ For this case the pseudo load vector becomes 


{F} k 


5[K] 

5V,_ 


{u} 


( 22 ) 


The analytical method consists of solving equation (20) for 9{u}/9v k using 
analytically computed derivatives for 9[K]/9V k (eq. (6)). Equation (20) is solved 
by the same solution algorithm used for solving equation (19), taking advantage of 
the fact that the factored form of [K] is available from the solution of equa- 
tion (19). Substitution of equation (6) into equation (20) leads to 


[K] 


5{u} 

5 \ 


a{f} 



NDDV NSDP 



a [k] 

9P . 

3 


9P. 5v. 

a vr^r {u} 

i k 


(23) 


Semi analytical Method 

The semianalytical method contains the analytical expression for the displace- 
ment derivatives (eq. (20) ) with finite-difference derivatives of the stiffness 
matrix (eq. (15) or eqs. (17) and (18)). If equations (15) and (20) are used, the 
following expression results: 


9{u} _ 9{f } 

9V, 9V_ 

k k 



(24) 


where 9{f}/9V k can be calculated analytically or by finite differences and the term 
in parentheses may be obtained by finite differences from equation (15) or by the LSK 
processor (eqs. (17) and (18)). Derivatives are calculated by creating individual 
pseudo load vectors {F } k and solving equation (24) using the factored [K] matrix 
from equation (19). This formulation is a significant simplification of the analyt- 
ical approach of reference 4 in that it avoids the need for element-dependent manipu- 
lations to handle bending- type elements. It also eliminates the complications in the 
analytical method when design-variable linking is included (eqs. (1) and (23)). 


Finite-Difference Method 

The simplest method to implement, but the most time consuming computationally 
(especially for large finite-element models), is the finite-difference method. In 
addition, the accuracy of the finite -difference method depends upon the perturbation 


One practical case when this term is not zero is when thermal loads are 
included. 
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step size. In this method, the original structure is analyzed; the structure is then 
modified by perturbing a design variable and reanalyzed* The displacement derivative 
is 


5{u} 

5v, 





(25) 


where {u}^ n is the displacement vector due to a perturbation in the kth independent 
design variable and {u}° is the displacement vector of the original structure* 


STRESS DERIVATIVES 
Analytical Method 

Element stresses {a} are related to joint displacements and element tempera- 
tures by the equation 

{a} = [S] {u} - [G]{a}(T - T ) (26) 

o 


where [S] is the stress-displacement matrix, [G] is the stress-temperature 
matrix, {a} is the vector of thermal expansion coefficients, T q is the stress- 
free element temperature, and T is the actual element temperature* Upon differ- 
entiation with respect to the kth independent design variable, the general expression 
for stress derivatives is 


a{a} 

av. 


NDDV NSDP 



5[S] 

5P. 

D 


dP. 5v. 

5v. 5V, 
i k 


(27) 


Equation (27) neglects any dependence of temperature on the design variables. 

For rods (E23), membranes (E31 and E41), and shear panels (E44), the matrix [S] 
is independent of the section-property design variables; thus, the second term on the 
right side of equation (27) is zero. For bending-type elements, beams (E21) and 
plates (E32, E33, E42, and E43), the stress-displacement matrices are dependent on 
the element cross-sectional geometry, and this dependence must be included in the 
stress derivative calculations* This was done in reference 4 for the channel-section 
beam in figure 2 by using specific, user-defined subroutines. 


Semi analytical Method 

In this method (reported in ref. 11), the perturbed displacement vector {u}^. n 
is approximated by a truncated Taylor series expansion as follows: 



(u}° 


5{u} 

av,_ 


(v v n - O 


(28) 


1 1 



where Stul/SV^ is calculated with equation (24). The perturbed stress vector is 
given by 


{a> k n = [S] k n (u} k n - [G]{a}[G]{a}(T - T q ) 
and the derivative is 


(29) 


5{cf} 

5V. 


io) k n - {of 



where 


{a}° = [S]{u}° - [G] {a} (T - T ) 

O 


(30) 


(31) 


Finite-Difference Method 

The stress derivatives are calculated in a manner analogous to the displacement 
derivatives. The structure is modified and reanalyzed and the derivatives are 
approximated as 


5{q} 

5v 


{oh n - (a}° 
k 

n o 

V - V, 
k k 


(32) 


where {a} k n is the stress vector due to a perturbation in the kth independent 
design variable calculated with equation (29) and {a}° is the original stress 
vector. 


DERIVATIVES OF VIBRATION AND BUCKLING EIGENVALUES AND EIGENVECTORS 

Methodologies for calculating derivatives of eigenvalues and eigenvectors 
have been presented in references 4, 6, 8, and 9. The methods which have been 
implemented in EAL are an analytical method (ref. 8) and the finite-difference 
method. 


Analytical Method 

This method is implemented in EAL using runstreams based on the development in 
reference 4 and the theory in reference 8. The matrix equation for free vibrations 
is 


1 2 


([K] - = 0 


(33) 



where 4}j and u. are, respectively, the vector of the mode shape and the scalar 
frequency corresponding to vibration mode j, [M] is the mass matrix, and [K] is 
the stiffness matrix* Normalizing the modes with respect to the mass matrix gives 


4) t [m] 4> = [i] 


(34) 


where [I] is the identity matrix. Differentiating equation (33) with respect to 
V k gives 



\ 


a{4>> . 


“j 


3 _ 


da). 


ay 


-[M] {<!>}. - 


a [k] 
ay 


4) . + to 


2 a [mi 


ay 


4}. 


(35) 


Premultiplying equation (35) by 


(4>} t 


and using equations (33) and (34) gives 



r., T a [K] <% 2 /aX T a[M] /Al 


K 


ay 


(36) 


Since [K] - (Oj^[M] is singular, a direct solution of equation (35) is not 
attempted. The value of one component of d{<J>}j/dV k is fixed, yielding a partic- 
ular solution {Q} j to equation (35). This is accomplished in EAL by identifying 
the component of the eigenvector with the largest absolute value and constraining to 
zero the corresponding component of the eigenvector derivative. The eigenvector 

is the complementary solution to equation (35). Thus, the expression for the 
eigenvector derivative is 


a {<{>> . 

“av = {Q} j + c{<})} j (37) 


The value of the multiplier C is obtained by substituting the expression for 
a{iJ)}j/5V k from equation (37) into the following: 




[M]- 


a{<t>}. 


ay 


i = -4). 


t a[M] 


av. 


4). 


(38) 


Equation (38) is the derivative of the egression for the normalization of modes with 
respect to the mass matrix (eq. (34)). The resultant equation for C is 



wnelj 


(39) 
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For buckling eigenvalue and eigenvector derivatives, the preceding equations 
apply directly, with the buckling mode shape {$}j and load X^ substituted for the 
vibration mode shape {<J)}j and frequency u) j 2 and with the negative of the geometric 
stiffness matrix -[Kg] substituted for the matrix [M] . In the present implemen- 
tation, derivatives of the matrices [K] , [M] , and [K g ] were calculated by finite 

differences. 


Finite-Difference Method 

The equations for finite-difference derivatives of eigenvalues and eigenvectors 
are analogous to those previously shown for displacement and stress derivatives 
(eqs. (25) and (32)). When the finite-difference method is used to obtain deriva- 
tives of eigenvectors, careful attention must be paid to accuracy of the eigen- 
vectors. Problems arise because eigenvalue routines generally base convergence on 
accuracy of the eigenvalues, whereas the associated eigenvectors may not be as 
accurate. Hence, when differencing the eigenvectors, further errors are introduced 
which may make the derivatives inaccurate. 


IMPLEMENTATION OF DERIVATIVE CAPABILITY IN EAL 

An outline of the implementation of the derivative capability in EAL is given in 
the appendix. For brevity, only the implementation of the semi analytical method is 
discussed in the appendix, which includes descriptions of the overall procedure, 
descriptions of the key portions of the EAL runstreams, and a listing of the input 
for the problem of the box beam. 


NUMERICAL EXAMPLES AND RESULTS 

Four example problems were investigated to verify the derivative implementation 
and to compare the methods. The problems were a swept wing, a box beam, a stiffened 
cylinder with a cutout, and a space radiometer-antenna truss. Comparisons of solu- 
tion times for the various derivative calculation methods are given, and a conver- 
gence study was performed for the wing problem to verify the accuracy of the 
finite-difference results and to establish an appropriate perturbation in the design 
variables for the remaining finite-difference calculations. Analytical derivatives 
(when available) were used to assess the accuracy of the other methods. All numer- 
ical calculations were performed on the CDC® CYBER 175 computer under NOS 1.4. 


Swept Wing 

Optimization of the swept wing shown in figure 3 has been investigated and 
reported by several researchers (refs. 17 to 19). This example was a moderately 
complex representation of a wing modeled by rods (E23), triangular membranes (E31), 
and shear panels (E44) . Since bending elements were not used, the calculation of 
derivatives with the analytical method was straightforward, and the results were used 
to check the accuracy of the semianaly tical and finite-difference methods. The 
geometry and node numbering are shown in figure 3; nodal coordinates are listed in 
table 3 and design variables are described in table 4. In this problem, 32 design 
variables controlled the section properties of 150 elements. The wing was aluminum, 
had a constrained root, and was subjected to two load conditions. (See table 5. ) 

Load condition 1 was approximately equivalent to a uniform pressure loading of 
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0.556 psi. Load condition 2 had the same total load, but the distribution was 
changed to move the center of pressure forward. The wing was symmetric with respect 
to the x-y plane (u and v displacements were zero on the plane z = 0) , and hence 
only half the wing was modeled. The model consisted of triangular membrane elements 
for the skin, rod elements for the spar caps, and shear panels for the rib and spar 
webs. (See fig. 4.) Forty-four transverse rod elements (not shown) were added at 
the vertical edges of the shear panels to provide the necessary stiffness throughout 
the depth of the wing. The thicknesses t i were 0.2 in. for E31 elements 1 to 24, 
0.1 in. for E31 elements 25 to 60, and 0.2 in. for all E44 elements. The cross- 
sectional area A i was 0.02 in 2 for E23 elements 1 to 20 and 0.2 in 2 for E23 
elements 21 to 64. 

A study was performed to determine the perturbation step size needed to obtain 
sufficiently accurate results for the finite-difference and semianalytical methods. 
The results are shown in table 6(a). The value of the largest displacement deriva- 
tive value (5 w/ 5V 14 at joint number 41 for load condition 1) was used to judge 
convergence. As indicated in table 6(a), the finite-difference result, with a 
1 -percent increment in the design variable, was within 1 percent of the analytical 
results. Use of perturbations smaller than 0.01 percent led to degradation of 
accuracy. The semianalytical method was capable of duplicating the analytical result 
with a 1 -percent perturbation in the design variable. Although the acceptable per- 
turbation sizes are generally problem dependent, this dependence was not signifi- 
cant. In all subsequent calculations, perturbations of 1 percent were used. 

Comparisons of solution times for displacement derivatives are also presented in 
table 6(b). The finite— difference method was slowest (700 sec), the semianalytical 
method required 419 sec without the LSK processor and 161 sec with the LSK processor, 
and the analytical method required 135 sec. Hence, the use of the LSK processor 
makes the semianalytical method significantly more effective and makes it competitive 
with the analytical method. 


Box Beam 

The next problem was a cantilevered box beam (fig. 5). This example was used to 
verify displacement derivative runstreams for all the EAL structural elements listed 
in table 2. 

The geometry and joint numbering are shown in figure 5, and elements are illus- 
trated in figure 6. The thickness of all two-dimensional elements was 0.1 in. The 
cross-sectional area of the rod elements was 1.0 in 2 , and the beam elements had tube 
sections with an inner radius of 2.0 in. and an outer radius of 2.5 in. The mate- 
rial was aluminum with the following material properties: p = 0.096 lb/in 3 , 

E = 10.6 x 10 psi, and v = 0.3* The applied load was composed of 10 000-lb forces 
in the positive z direction at nodes 27 and 28. The cross-sectional area or thick- 
ness of each group of elements was considered to be a separate design variable. The 
problem had 10 dependent design variables and 9 independent design variables. The 
inner and outer radii of the tube elements were linked to an independent variable r 
by the equation 



(40) 
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Displacement derivatives with respect to each independent design variable were 
calculated with the finite-difference method and the semianalytical method. 2 Results 
shown in table 7(a) indicate run times of 133 sec for the finite-difference method 
and 63 sec for the semianalytical method with the LSK processor. The reduction in 
solution time for the semianalytical method with LSK over the finite-difference 
method was smaller for this example than for the swept wing because the swept wing 
problem had more independent design variables (32 instead of 9). Values shown in 
table 7(b) reveal which design variables would be most effective to increment to 
reduce tip deflections. As expected, the upper and lower covers (coupled plate 
elements E43 and E33) had the greatest effect on limiting tip deflections. In 
addition, since the E32 and E42 elements had only bending stiffness and the bending 
loads in the covers were negligible, the derivatives with respect to thickness of 
these elements were negligible. The results in table 7(b) suggest that derivatives 
can be used directly to guide improvements in a structural design apart from their 
use in formal optimization. In some instances substantial improvements may be 
realized (ref. 20). 


Stiffened Cylinder With Cutout 

Optimization studies of a cantilevered stiffened cylinder with a cutout (fig. 7) 
were reported in references 4 and 21. Node point locations are given in table 8. 
Element types and locations are shown in figure 7(b). The cylinder model had 
352 degrees of freedom and was stiffened by 5 equally spaced rings along the length 
and 16 equally spaced stringers around the circumference. Rotations about the Y-axis 
at all node points and all translations at node points 1 to 1 6 were constrained to 
zero. The rings and stringers were modeled by beam (E21 channel sections) and rod 
elements (E23). Rectangular panels between rings and stringers were modeled by 
membrane elements (E41). The cross-sectional area of all rod elements was 0.646 in 2 , 
the thickness of the rectangular membranes was 0.0394 in., and the channel- 
section dimensions are shown in figure 7(b). The material was aluminum alloy with 
E = 10.8 x io 6 psi and v = 0.3. The independent design variables for this problem 
included the area of the rods, the thickness of the membrane elements, and a scale 
factor which controlled the cross-sectional dimensions of the beam elements. The 
loading consisted of two equal and opposite concentrated axial forces of 20 000 lb 
applied at the free end. 

Derivatives of displacements, stresses, vibration frequencies and mode shapes, 
and buckling loads and mode shapes were calculated. Computer times for calculating 
derivatives are presented in table 9. Use of the semianalytical method without the 
LSK processor for displacement derivatives resulted in a reduction in solution time 
of about 25 percent over the finite-difference method. Using the semianalytical 
method with the LSK processor resulted in a reduction in solution time of 28 percent 
over the finite-difference method. The reason solution times for the semianalytical 
methods were so similar for this case was that the ratio of the number of independent 
to the number of dependent design variables was only 3 to 5. As this ratio 
decreases, the advantage of the LSK processor increases. (In the problem of the 
swept wing, the ratio was 32 to 150 and the LSK processor had a much larger effect.) 

Stress derivatives were calculated by finite differences (122 sec) and semi- 
analytically (95 sec). The savings in solution time for finite differences relative 


2 


Derivatives for the membrane and rod elements were computed with the analytical 


method. 
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to the semi analytical method were not as dramatic in this problem as in the problems 
of the swept wing or the box beam because of the smaller number of design variables. 

Computer times for calculating derivatives of vibration frequencies, vibration 
mode shapes, buckling loads, and buckling mode shapes with the finite-difference and 
analytical methods are also shown in table 9. For the vibration problem, the 
analytical method was 57 percent faster than finite differences, and for the buckling 
problem it was 37 percent faster. 


Antenna Truss 


A finite-element model of a 180-ft radiometer-antenna reflector is shown in 
figure 8. The reflector, described in reference 22, was made up of tetrahedral truss 
modules, and the model consisted of 109 structural node points (table 10) and 420 rod 
elements. The structure was composed of graphite-epoxy composite with an effective 
modulus of elasticity of 10.6 x i psi and a coefficient of thermal expansion of 
0.13 x icT 6 per degree Fahrenheit. The antenna was subjected to thermal loading 
corresponding to Earth orbit at an altitude of 216 miles in an Earth-facing orien- 
tation. The thermal loading consisted of a combination of solar. Earth -ref lected 
(albedo), and Earth-emitted heat flux. A transient thermal analysis of the structure 
for a complete orbit was performed (ref. 23) and node point temperature differences 
from a worst-case condition (largest temperature gradients, table 11) were used for 
the present calculations. Three design variables were used for this problem. They 
were the cross-sectional area of the elements in the upper surface (0.2530 in 2 ), the 
area of the diagonal elements joining the upper and lower surfaces (0.1741 in 2 ), and 
the area of the elements in lower surface (0.2530 in 2 ). Both stress and displacement 
derivatives were calculated, and for this problem the applied load vector {f} was a 
function of the design variables. (See eq. (20).) 


Solution times for displacement and stress derivatives are summarized in 
table 12(a). The semianaly tical method was 39 percent faster than the finite- 
difference method for displacement derivatives and 29 percent faster for stress 
derivatives for this problem. As indicated by the tabulated derivatives in 
table 12(b), the most effective way to reduce the center deflection of the reflector 
would be to decrease the areas of the upper-surface elements. 


CONCLUDING REMARKS 

A capability for computing structural-sensitivity derivatives has been imple- 
mented in the general purpose, finite-element computer program denoted the Engineer- 
ing Analysis Language (EAL) System. This paper presented the development of 
runstreams which calculate derivatives of displacements, stresses, vibration fre- 
quencies and mode shapes, and buckling loads and mode shapes with respect to 
structural variables which include thicknesses, areas, and moments of inertia. 

Linear and nonlinear design-variable linking representations were included. Three 
methods for computing the derivatives were documented: analytical (which calculates 

exact structural-sensitivity derivatives from the governing equations), 
semianalytical (which combines finite-difference derivatives of mass, stiffness, and 
geometric stiffness matrices with analytical expressions for derivatives), and finite 
differences. The derivative capability was demonstrated for the following four 
structures: a swept wing, a box beam, a stiffened cylinder with a cutout, and a 

space radiometer-antenna truss. Comparisons of solution times for the various 
methods were described. A convergence study was performed for the swept wing to 
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verify the accuracy of the finite-difference results and to establish an appropriate 
increment in the design variables for use in the subsequent finite-difference calcu- 
lations* Analytical derivatives of displacements for the wing and of vibration fre- 
quencies of the cylinder served as benchmarks to assess the accuracy and efficiency 
of the finite-difference and semi analytical methods* Results indicated that the 
semianalytical method for calculating displacement and stress derivatives was 
efficient, general, and straightforward to implement* Further, the semianalytical 
method was not adversely affected by the presence of bending elements or design- 
variable linking* 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 8, 1984 
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APPENDIX 


SYSTEM OF EAL RUNSTREAMS FOR CALCULATING STRUCTURAL-SENSITIVITY 
DERIVATIVES WITH SEMI ANALYTICAL METHOD 

Overall Procedure 

This appendix describes a system of EAL runstreams used to calculate structural- 
sensitivity derivatives with the semianalytical method. A flowchart for the overall 
system is shown in figure 9. Runstreams or groups of runstreams presented at the end 
of this appendix are keyed by numerals to specific portions of flowcharts shown in 
figures 10 to 12. The semianalytical method uses equations (15) and (16) to calcu- 
late derivatives of the stiffness and mass matrices and uses equations (24) and (28) 
to (30) to calculate the displacement and stress derivatives. 

As shown in figure 10, derivatives of the stiffness, mass, and geometric stiff- 
ness matrices are calculated by finite differences. The procedure is as follows: 
loop over the independent design variables, perturb a design variable, link design 
variables to form vector of dependent design variables, loop over the dependent 
design variables and modify the structure, calculate new stiffness and mass matrices 
and geometric stiffness matrix (if appropriate), differentiate the resulting matri- 
ces, and restore the original value of independent design variable V fc . The numerals 
beside specific boxes correspond to runstreams described subsequently. 

Displacement derivatives (fig. 11) are calculated by looping over the indepen- 
dent design variables, building a pseudo load matrix (eq. (20) ) and solving equa- 
tion (24). The acquired displacement derivatives are then used to approximate a 
perturbed displacement vector {u}^ n ( e< T* (2 8) ) and then to approximate a perturbed 

stress vector (eq. (29)). The stress derivatives are then calculated by finite 
differences (eq. (30)). 

The flowchart which illustrates the calculation of derivatives of vibration and 
buckling eigenvalues is shown in figure 1 2. The outer loop is on the mode for which 
derivatives are required. The two inner loops are both over the independent design 
variables. The first inner loop calculates the derivatives of the vibration fre- 
quencies and buckling loads with equation (36) and then builds a pseudo load vector 
(right-hand side of eq. (35)). A particular solution of equation (35) is calculated 
and used in the second inner loop to calculate derivatives of the eigenvectors with 
equations (37) and (39). 

Descriptions of the key runstreams used in the system of figures 9 to 1 2 
(particularly those which are problem dependent) are given in the following table: 
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DESCRIPTION OF KEY RUNSTREAMS 


Runstream name 


Identification no, 
in figs. 9 to 1 2 


Description 


DRIV GRD4 


I NIT MODL 


SET PARA 


Driver . runstream which calculates derivatives semianalyti- 
cally? the entire system can be run by the statement 
♦DCALL (DRIV GRD4) 

Sequence of EAL processors (TAB, ELD, TAN, and AUS) which 
describes the initial structural model and load cases? the 
example shown in the runstream listing is for the box beam 

Initializes following parameters used by system to determine 
specifics of derivative information required? 


NLST 

Number of load sets 

MODT 

Number of modes for which derivatives are 
computed by method of reference 8 

MODE 

Same as MODT 

OLIB 

Output library for intermediate buckling and 
vibration derivatives 

OTIB 

Output library for intermediate stress or 
displacement derivatives 

NPOL 

Degree of polynomial used in linking design 
variables (eq. (1)) 

DISP 

If equal to 1, displacement derivatives are 
calculated 

STRD 

If equal to 1, stress derivatives are 
calculated 

BUCD 

If equal to 1, buckling load and mode shape 


derivatives are calculated; if equal to 2, 
only buckling load derivatives are 
calculated 

VIBD Same as BUCD for vibration derivative 

NJ Number of joints in model 

NDV Number of independent design variables 

NDVD Number of dependent design variables 

NDVA Number of rod and membrane dependent design 

variables (E23, E31, E41, E44) 

NBDV Number of beam dependent design variables 

(E21) 

NPDV Number of plate dependent design variables 

(E32, E33, E42, E43) 

NDVX Number of shape-dependent design variables 

DR Incrementing factor for design variables for 

finite-difference calculations 

DE21 , DE23, l Number of dependent design variables for each 

..., DE44 J element 

NE21, NE23,[ Number of independent design variables for 

..., NE44 r each element 

MODE NUM Name of table of mode numbers for which 

TABLE derivatives are taken 
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Runstream name 


Identification no. 
in figs . 9 to 1 2 


Description 


4 Set up tables denoted DVA DFN, DVP DFN, and DVB DFN which define 

membrane, plate, and beam dependent design variables, respec- 
tively; element types are defined by a numeral: E23-1, E31-2, 

E32-3, E33-4, E41-5, E42-6, E43-7, E44-8, and E21-9? the type 
of beam used in EAL in section 3.1.9 of reference 12 is denoted 
by a code as follows: TUBE - 1, BOX - 2, TEE - 3, ANG - 4, 

WFL - 5, CHN - 6, ZEE - 7, GIVN - 8, DSY - 9; beam design 
variable numbers are designated by the following codes: 


b 1 - 1, t 1 - 2, b 2 - 3, t 2 - 4, b 3 - 5, t, - 6, I - 7, 

«1 - 8, I 2 - 9, a 2 - 10, a - 11, f - 12, f 1 - 13, - 14, 

Z 2 “ 18, 0 - 16, q 1 - 17, q 2 - 18, q^ - 19; each of the 

variables shown is described in section 3,1.9 of reference 12; 

the order in which the dependent design variables are listed 
(starting with the membrane, then plate, and then beam design 
variables) are the actual locations of each design variable in 
the dependent design variable vector (see INTL DESV) ; the only 
restriction is that the structural parameters of beams with the 
same section property number must be listed together in the 
DVB DFN table; this permits the linking of several section 
properties of a beam; this could be used to facilitate the 
scaling of an entire beam cross section by one design variable 
or scale factor 


INTL DESV 


Initializes design variables by forming a table of initial design 
variables called DESV CNMN whose only limitation is that when 
linked it will produce the same structural parameters as the 
initial structural model (see INIT MODL) 


COEF LINK 


Series of tables which define the coefficients for nonlinear 

design-variable linking as expressed by equation (1) and rewrit- 
ten as {v} = {Cj + [C^Hv} + [C ] { V } , ..., + (C Hv* 1 } 
where {C Q } is represented as table COEF DV 1 0, m (C 1 ] is 
COEF DV 1 1, [C 2 1 is COEF DV 1 2, ..., [C M ] is COEF DV 1 M; 

each block of the above tables represents a column in the 
[C] matrices; this runstream must be included in the system of 
runstreams even if linking is not used 
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Runstream Listings 

This section of the appendix contains a listing of the EAL runstreams used to 
calculate derivatives of displacement, stress, and vibration and buckling eigenvalues 
and eigenvectors with the semianalytical method. 


♦ char $*i 

(♦ECHO 1*2*3, 4 
♦ABORT 1 
♦XQT U1 
$ 

( 1 

S 

♦ ( DR I V GRD4) ENDDG4 

(DRIVER PROGRAM TO CALCULATE DU/DCAPV BY THE S EM I -AN AL YT IC AL METHOD 
♦DC ALL ( 1NIT MODL) 

♦XQT DCU 
DUPLICATE 1 2 
DUPLICATE 1 3 
DUPLICATE 1 4 
DUPLICATE 1 5 

♦ DCALL ( S ET PARA) 

♦ DC ALL ( D 6S V DFN ) 

♦DCALL ( INTL DESV) 

♦ DCALL ( COE F LINK) 

♦DCALL (CALC GRD4) 

♦DCALL ( PRT SUB) 


♦RETURN 

♦ ENDDG4 

( 

2 

$ 

♦ ( I N 1 T MODL) ENDINM 

$$ (( INITIAL MODEL STRUCTURAL DEFINITION 

♦XQT TAB 
START 2B 


TITLE 1 BOX BEAM GRADIENT CALCULATION 
TEXT 

•GRADIENT CALCULATION OF A BOX-BEAM WHICH USES 
•E21*E23*E31*E32*E33*E41*E42*E43*AND E44 ELEMENTS 
(INPUT JOINT LOCATIONS 
JLOC 

1 0 . 0 . 0 . 0 . 10 . 0 . 217 

2 60. 0. 0. 60. 10. 0. 

15 0. 0. 10. 0. 10. 10. 217 

2 60. 0. 10. 60. 10. 10. 

(SPECIFY MATERIAL PROPERTIES 
MATERIAL CONSTANTS 
1 10.6E+6 .3 .096 

(DESCRIBE MATERIAL REFERENCE FRAME 
MREF 

FORMAT-2 

1 * 1 * 1000 . * 0 . * 1000 . 

2 * 1 * 1000 . * 10 . * 1000 . 

(BEAM ELEMENT SECTION PROPERTIES 
BA 

TUBE 1* .25* .3 
(ROD ELEMENT AREAS 
dC 

1 * 1 . 

(MEMBRANE ANS PLATE ELEMENT THICKNESSES 
SA 

NMAT-1 

1* .1$E31 ELEMENTS 
2* . 1 (E 32 ELEMENTS 
3* • 1 (E 33 ELEMENTS 
4* .1$E41 ELEMENTS 
5* • 1 $E42 ELEMENTS 
6* .1(643 ELEMENTS 
(SHEAR PANEL THICKNESSES 
SB 

1 * .1 

(CONSTRAINT DEFINITION 
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CON* 1 

ZERO 1,2>3#4,5j6J1,2 
ZERO 1,2>4,5, 6)15/16 
♦XQT ELD 

SBEAM ELEMENT DEFINITION 
E21 

NMAT *1 
NSECT-1 
NRE F * 1 
13 16 

15 17 1 6 
NREF *2 

2 4 16 

16 18 1 6 

SRQD ELEMENT DEFINITION 
E23 

NMAT *1 
NSECT-1 

1 15 1 1 7 2 

2 16 1 1 7 2 

$TRI ANGULAR MEMBRANE ELEMENT DEFINITION 
E31 

NMAT »1 
NSECT-1 

8 10 7 

7 9 10 

9 11 10 

10 12 11 
12 14 11 

11 13 14 

SUNCOUPLEO TRIANGULAR PLATE ELEMENT DEFINITION 
E32 

NMAT *1 
NS ECT *2 

8 10 7 
7 9 10 

9 11 10 

10 12 11 

12 14 11 
11 13 14 

SCQUPlED TRIANGULAR PLATE ELEMENT DEFINITION 
E33 

NMAT* 1 
NS EC T* 3 

13 2 

2 4 3 

4 6 3 

3 5 6 

5 7 6 

6 8 7 

$M QUADRILATERAL MEMBRANE ELEMENT DEFINITION 
E41 

NMAT -1 
NS ECT ■ 4 

21 23 24 22 131 

SUNCOUPLEO QUADRILATERAL PLATE ELEMENT DEFINITION 
E42 

NMAT -1 
NS ECT "5 

21 23 24 22 131 

SCOUPLEO QUADRILATERAL PLATE ELEMENT DEFINITION 
E43 

NMAT *1 
NS EC T *6 

15 17 18 16 131 

SSHEAR PANEL DEFINITION 
E 44 

NMAT *1 
NSECT-1 

1 3 17 15 1 6 1 

2 4 18 16 1 6 1 
♦ XQT TAN 

SAPP L I E 0 LOADS DEFINITION 
♦XQT AUS 

SYS VEC } APPLIED FORCE 1 

1*3 J J»27# 28 ; 10000, 10000. 

♦RETURN 
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* ENDINM 

$ 

S 3 

S 

♦(SET PARA ) ENDPAR 

SS S$ SET CONTROL PARAMETERS 


INLST-1SNUMBER OF LOAD SETS 
ILCAS-ISNUMBER OF LOAD CASES 
IM JDT »2$T JTAL NUMBER OF MODES 

1DLIB-15SOUTPUT LIBRARY FOR INTERMEDIATE DYNAMIC DERIVATIVES 
IOTIB-15SOUTPUT LIBRARY FOR INTERMEDIATE STRESS OR DISP. DERIVATIVES 
IULIB-15 
l MOD E* 2 

INPOL-1SDEGREE OF POLYNOMIAL LINKING 
IDI5P-1SDIS PLACEMENT D ERI VER AT IV ES 
ISTRD" 1SSTRESS DERIVATIVES 
IBUCD-1S6UCK DERIVS.(LOADEFREO) 

! V IBD* 1$VI 8, DERIVS. (LOADEFREQ) 

1NJ-26SN0. OF JOINTS 

I NDV ■9$N0. OF INDEPENDENT DESIGN VARIABLES 
INUDV-NDV 

SNDVX-OSNO.OF DEPENDENT SHAPE VARIABLES TO BE LINKED 

INDVA-A $NO. OF DEPENDENT MEMBRANE -TYPE DESIGN VARIABLES TO BE LINKED 

!NLNK*10S NO, OF LINKED DESIGN VARIABLES 

!NDV1-"NDV" 

INBDV-2S NO. OF B t AM DEPENDENT OESIGN VARIABLES TO BE LINKED 
lNPDv-AS NO. OF PLATE DEPENDENT DESIGN VARIABLES TO BE LINKED 
1NDVD-NDVX+NDVA+NBDV+NPDV 
INDVI-NDV 

i DR* . 001$ FACTOR FOR INCREMENTING DESIGN VARIABLE 

IDE 23* 1SNQ • OF E23 D.V.'S 

1DE31-1S NO. OF E3I D.V.'S 

IDE32-1S NO. OF E32 D.V.'S 

IDE 33 s IS NO, OF E33 D.V.'S 

IDEA1-1S NO. OF EA1 D.V.'S 

l DE A2 ■ 1 $ NO. OF EA2 D.V.'S 

l DE A 3* 1 $ NO. OF E A 3 D.V.'S 

IDE AA*1$ NO. OF EA A D.V.'S 

IDE21-2S NO. OF E21 D.V.'S 

I NE 2 1* 1 

IINX1*DS,A,1,1(1# JDF1, BTAB,1,8) 

IINX2-DS, 5,1, 1(1, JDF1,BTAB, 1,8) 

!INX3-D5,6,1,1(1, JDF1,BTAB,1,6> 

IINXA-DS,7,1,1(1,JDF1,BTAB,1,8) 

!INX5-DS,8#1,1(1, JDF1, BTAB,1,8) 

IINX6-DS,9,1,1(1,JDF1,BTAB,1,8) 

$ 

STABLE OF MODE NUMBERS 

$ 

♦XQT AUS 

TABLE (RMODE* 2, TYPE-0, NI»1,NJ -"MODE"); MODE NUM 

J-lJl 

J-2;2 

$ 

STABLE OF UNCONSTRAINED DEGREES OF FREEDOM 

S 

TABLE(RM0DE-2,TYPE-0,NI-1,NJ-6)J IN EX 
J-1;"INX1" 

J-2;"INX2*' 

J-3;"INX3" 

j-a;"INXA” 

J«5)"INX5" 

J-6;"1NX6" 

IINX1-FREE ( ) 

I1NX2-FREEU 
II NX3*FREE ( ) 

IINXA-FREEU 
UNX5-FREEO 
1 I NX6- FREE( ) 


♦RETURN 

♦ ENDPAR 
S 

S A 

S 

♦ ( DE S V DFN ) ENDDVD 

SS S$ DEPENDENT OESIGN VARIABLE DEFINITION 

♦XQT AUS 
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SMEMBRANE-TYPE DESIGN VARIABLES 

TABLE (RM0DE-2»TYPE«0#NI-4,NJ -"NOVA" JJDVA DFN 
J*l» M NCVA"$ELEflENT TYPE GROUP NO. NO. OF ELEN. SECT. NO. 

1 1 12 1 

2 16 1 

5 1 3 A 

6 1 12 1 

SPLATE-TYPE DESIGN VARIABLES 

TABLE (RMODE-2,TYPE»0,NI»4,NJ» M NPDV">JDVP DFN 
J«1>"NPDV" SELENE NT TYPE GROUP NO. NO. OF ELEN. SECT. NO. 

3 1 6 2 

A 1 6 3 

6 1 3 5 

7 1 3 6 

S8EAM-TYPE DESIGN VARIABLES 

TABLE (RNDDE«2»TYPE»0>NI»5»NJ«"NBDV”);DVB OFN 
J*1»"NBDV"$ELEM. TYPE GROUP NO. NO. OF ELM. SECT. NO. D.V. NO. 


1 1 24 11 

1 1 24 12 

♦RETURN 

♦ ENDDVD 

S 

$ 5 

$ 

♦tINTL DESV) ENDINO 

Si SS INITIALIZE DESIGN VARIABLES 

♦XQT AUS 


TA3LE <NJ«"NDV'') ) DESV CNMNSDATA SET CONTAINING INITIAL VALUE OF DESIGN VARIABLES 

J * 1 J 1 • 

J*2> 6; • l 
J*9; .25 
♦RETURN 


♦ ENDIND 

$ 

$ 

$ 

♦( CQEF LINK) ENOCLK 

SS $$ MATRIX OF COEFFICIENTS FOR LINKING 


♦XQT AUS 

TABLE(NJ-”NDVO») JCOEF DV 1 1 
BLOCK 1 
i-i; j-isi. 

BLOCK 2 
I*i;j*2;i. 

BLOCK 3 
i«i; 

BLOCK 4 
i-li 
BLOCK 5 
i*i; 

BLOCK 6 
1*1 ; J*6; l . 

BLOCK 7 
i-l; J«7;i. 

BLOCK 8 
1*1; J*0 jl. 

BLOCK 9 
l*i; J-9;i. 

1*1; j-iojl. 

TABLE(NJ*"NDVD") JCOEF DV 1 0 


J* If 9) 0. 

J* 10 J • 05 
♦RETURN 

♦ ENDCLK 

S 

, 7 
♦(CALC GRD4 ) ENDGD4 


ss $s CALCULATE DERIVATIVES OF U W.R.T. DESIGN VARIABLES USING THE SEMI-ANALYTICAL METHOD 

♦DC ALL (CALC KM) 

♦ DC ALL ( DKDV DMDV ) 

♦ JNZ ( D ISP> 262 ) 

♦JNZ ( STRD# 262 ) 

♦JUMP 265 
♦LABEL 262 

♦DC ALL ( CALC DUDV) 

♦ LABEL 26<r 

♦JZ (STRD*265)$TEST TO SEE IF STRESS DERIVATIVES ARE REQUIRED 


25 



APPENDIX 


♦DCALL ( 5TRS DERV ) SCALCUL ATE STRESS DERIVATIVES 
♦LABEL 265 

$TEST TO SEE IF BUCKLING DERIVATIVES ARE REQUIRED 
♦JZ ( BUCD#266 ) 

1N-BUCD 

* JNZ > -l (N > 267 ) 

INAM9- •BUCK 
INAM7- 1 BKVL 
!NAM8* 1 3KMD 
♦DCALL (NELS METH) 

♦JUMP 266 
♦LABEL 267 
♦JNZ>-1<Nj269) 

INAM9- 'BUCK 
INAM7" 'BKVL 
1NAM8-'BKMD 
♦DCALL (NELS METZ) 

♦LABEL 269 
♦JNZj-1 (N* 273 ) 

IN AM9* ' BUCK 
INAM7- ' BKVL 
iNAMB-'BKMD 
♦DCALL (MODL METH) 

♦JUMP 266 
♦LABEL 273 
♦JNZ, -1<N, 285) 

INAM9- ' BUCK 
1NAM7- ' 8KVL 
1 NAM8- * BKMD 
♦DCALL (NELS METH) 

♦JUMP 266 
♦LABEL 285 
♦XQT EXIT 
♦LABEL 266 

♦JZ(VIBD,268)$TEST TO SEE IF VIBRATION DERIVATIVES ARE REQUIRED 
tN-V IBD 

♦JNZ,-1(N*270) 

INAM9- 'VIBR 
iNAM7- * VBVL 
1NAM8»VBMD 
♦DCALL (NELS METH) 

♦JUMP 268 
♦LABEL 270 
♦JNZ,-1(N>271) 

JNAM9- ' VIBR 
IN AM7- • VBVL 
!NAM8- ' VBMD 
♦DCALL (NELS MET2 ) 

♦JUMP 268 
♦LABEL 271 
♦ J NZ, -1 { N, 274 ) 

{NAM9- 'VIBR 
!NAM7- ' VBVL 
iNAM8- ' VBMD 
♦DCALL ( MOOL METH) 

♦JUMP 264 
♦LABEL 274 
♦JNZ,-1 (N, 275 ) 

IN AM 9- 'VIBR 
!NAM7- ' VBVL 
JNAM8- ' VBMD 
♦DCALL ( MUDL METH) 

♦JUMP 268 
♦LABEL 275 
♦XQT EXIT 
♦LABEL 268 


♦RETURN 

♦ ENDGD4 
$ 

$ 8 

$ 

♦(CALC KM) ENDCKMSPERFORM STATIC AND/OR DYNAMIC ANALYSES OF INITIAL STRUCTURE 

♦ XQT E 
♦XQT EKS 

♦ XQT K 


♦ JNZ t D I S P # 156) 
♦JNZ (STRD, 156) 
♦JUnF 157 
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♦LABEL 156 
♦XQT RSI 
! A*NLST 
♦LABEL 5666 
♦JLZ>— 1(A*5667 ) 

INNST-NLST-A 

♦XQT SS‘JL$CALCULATE DISPLACEMENTS 
RESET SET- M NNST” 

♦JUMP 5666 
♦LABEL 5667 
1A-FREEU 
INNNST»F RE E ( ) 

♦LABEL 157 
♦XQT U1 

INN3«T0C#N3(1,K, SPAR, MASK, MASK) 

INN W0C,N4(1,K, SPAR, MASK, MASK) 

♦JZ(STRD,150> 

IA-NLST 
♦LABEL 7666 
♦JLZ,-1(A,7667> 

1NNST-NLST-A 
INAnl-»STAT 
I NAM2" 1 DIS P 
INAM3*"NNST ,# 

InAMA-I 
! CiL I B * 1 

IN11— 1000 $SET Nl— 1000 FOR UNPERTURBED STRESS CALCULATIONS 

♦DC ALL (CALC S TRS ) SC ALCUL ATE STRESSES 

♦JUMP 7666 

♦LABEL 7667 

!A-FkEE< ) 

INNST *FREE U 

INAM1-FREEU 

1NAM2-FREEO 

INAM3-FREEU 

INAM^-FREEO 

I0LIB-0TI3 

♦LABEL 150 

♦ JNZ (DISP# 193) 

♦JNZ(STRD,193) 

♦XQT RSI 

♦XQT SSOL 
♦LABEL 193 

♦ J Z ( BUCO, 151 ) 

♦XQT GSF 
RESET EMBED-1 
♦XQT KG 

♦XQT EIGSPERFORM BUCKLING ANALYSIS 
RESET CONV * 1 . - 10 
RESET PROB-BUCK 
RESET INIT-"MMDB'»,NREQ-"MRQB" 

RESET NDYN-20 
♦XQT EIG 

RESET C0NV-1,-10,PR0B«8UCK,NREQ-"MRQB",NDYN«20 
♦LABEL 151 
♦J Z ( V I BD, 152) 

♦XQT E4SPERF0RM VIBRATION ANALYSIS 
RESET NMODES-"MMOD",NREQ«"MREG M 
♦LABEL 152 
♦XQT DCU 
♦JNZ(DISP,159) 

♦ JNZ ( S TR D, 159 ) 

♦JUMP 156 
♦LABEL 159 
IA-NLST 
♦LABEL 6666 
♦JLZ,-1(A,6667) 

INNS T»NLST-A 

SRENAME ALL DATA SETS AS INDICATED BELOW 

CHANGE 1 STAT DISP "NNST" 1, PREV DISP »NNST" 1 

♦JUMP 6666 

♦LABEL 6667 

♦LABEL 159 

!A-FREE( ) 

I NNST*FR EE ( ) 

♦LABEL 158 

CHANGE 1 OEM DlAG 0 0,DEMP D I AG 0 0 
CHANGE 1 K SPAR "NN3” "NN4 M ,KP SPAP M NN 3 ” 
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CHANGE 1 INV K 1 0,INV KP 1 0 
CHANGE 1 X IN V K 1 0,XINV KP 1 0 
*J 2 ( BUC D, 154 ) 

1NNN3*TQC,N3(1, KG, SPAR, MASK, MASK ) 

!NNNW0C,NM1,KG, SPAR, MASK, MASK) 

CHANGE 1 KG SPAR "NNN3" »NNNV»,KGP SPAR "NNN3" "NNN4» 

CHANGE 1 BUCK EVAL 1 1,PREV BKVL 1 1 

CHANGE 1 BUCK MODE 1 1#PREV BKMD 1 1 

♦LABEL 154 

♦JZ< VIBD,155) 

CHANGE 1 VIBR EVAL 1 1#PREV VBVL 1 1 

CHANGE 1 V IBR MODE 1 1#PREV VBMD 1 1 


♦LABEL 155 
♦RETURN 

♦ ENDCKM 
S 

% 9 

S 

♦ ( DKD V DM0 V ) ENODKHSCALCULATE DK/DV#Dii/D V# AND DKG/DV 

1IIV-0 

IIIEVNDV 
♦LABEL 803 


♦ JLZ,-1(HEV,804)SL0UP OVER DESIGN VARIABLES 
1 1 IV * 1 1 V + 1 

!DVP-DS,1,"IIV",1(1,DESV,CNMN,MASK,MASK) 

IDV-DR+DVP+DVP 
1AQ-1. /DR 
1AQ-AQ/QVP 
1 8Q*-1 •♦A 3 
♦XQT AUS 

TABLE, U(NJ-"NDV")JDESV CNMNS INCREMENT DESIGN VARIABLE 
OPE R“XSUM 
J« M H V M # "OV" 

♦DC ALL ( TRNS DESV) 

♦DC ALL (LINK POLY) 

♦XQT DCU 
PRINT 1 LNK DV 
♦XQT TAB 
UPDATED 

♦ J Z ( ND V A# 8 00 ) 

!IE V»NDVA 

♦ DC ALL ( OVA F DD ) $ UPDATES ALL MEMBRANE PROPS. TO CURRENT VALUE 
♦LABEL 800 

♦JZ(NPDV#801) 

IIEV-NPDV 

♦DC ALL { DVP FDD)* UPDATES ALL PLATE PROPS. TO CURRENT VALUES 
♦LABEL 801 

♦ J Z ( N B D V , 6 02 ) 

IIEVNBDV 

♦ DC ALL ( DVB FDD ) 1 UPDATES ALL BEAM PROPS. TO CURRENT VALUES 

♦label 802 

JNN3«T0C#N3t 1#KP# SPAR, MASK #M ASK) 

♦ XQT E 
♦XQT EKS 

♦ XQT K 

RESET OUTLIB- "OLIB" 

♦XQT AUS 

DEFINE KN- "0L16" K SPAR 
DEFINE K0*1 KP SPAR 
DK-SUM(»AQ" KN,"BG” KO) 

DEFINE A-DEMP DIAG 0 0 

DEFINE B-DEM DIAG 0 0 

DMDV DIAG 0 "I IV »-SUM< »AQ" B , M B Q M A) 

♦XQT DCU 

SPRINT 15 K SPAR 
SPRINT 1 KP SPAR 

CHANGE 1 DK MASK MASK MASK, DK DV SPAR »NN3» "IIV" 

SPRINT 1 DKDV SPAR "NN3" "IIV” 

♦JZ(BUC0,294) 

♦XQT RSI 

RESET NUK* H 0LI8 M 

♦ XQT S SOL 

RESET KLIB-»ULIB» 

RESET KILIB-I 
♦XQT GSF 
RESET EMBED-1 
♦LABEL 292 
♦ JNZ ( BUCD, 293 ) 
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♦JNZ< V1BD,293) 

♦JUMP 294 
♦LABEL 293 
♦XGT KG 
♦XQT AUS 

DEFINE KGN-KG SPAR ,, NNN3 M "NNN4" 

DEFINE KGO«KGP SPAR W NNN3” M NNN4 H 
DKG* S UM ( #, AQ" KGN, #, BQ M KGO ) 

♦XQT DCU 

CHANGE I OKG SPAR "NNN3" "NNNA'SOKG SPAR "NNN3" "IIV" 

♦LABEL 294 
♦XQT AUS 

TABLE, U(NJ* M NDV")JDE$V CNMNS SET ORIGINAL VALUE JF DESIGN VARIABLE 
□PER "XSUM 

j* n l iv M ; "DVR" 

INIEV* IIEV 
♦ JNZ { N IE V, 805 ) 

♦DC ALL ( TPNS DESV) 

♦DC ALL ( LINK POLY) 

♦LABEL BOS 

♦ DC ALL (DCU EFIL) 

♦JUMP 803 
♦LABEL 604 

!0I 1«FRE£ ( ) 

IDI2-FREE( ) 

IDI3*FREE( ) 

!DI4-FREE( ) 

1DI5-FREE< ) 

! D 16 * F RE E ( ) 

ID 1 7-FREE € I 
1DI8*FREE ( ) 
t DI 9-FREE C ) 

IDIIO-FREEO 
I D I 1 1 ■ FREE ( ) 

!DI 1 13*FRE E ( ) 

JDI14-FREEO 
JDI15-FREEU 
IDI16-FREEO 
(DI17-FREEO 
JDI18-FREEO 
!DI19*FREE< ) 

IDI20-FREEO 
1DI21«FREE ( ) 

IDI 22*FREE ( ) 

I D 1 2 3*FREE ( ) 

!DI24-FREE( ) 

IDI25-FREE ( ) 

♦RETURN 

♦ ENDDKM 

S 

s 

s 

♦(CALC OUDV) ENDDUVSCALCULATE DU/DCAPV 

1A-NLST 

♦LABEL 5666 

♦ J LZ,-1 ( A, 5667 ) SLOOP OVER LOAD SETS 
iNNST-NLST-A 

1 1 1 v«o 

IIIEV-NDV 
♦LABEL 7777 

♦JLZ,-1(IIEV,7778)SL00P OVER DESIGN VARIABLES 

iiiv-imi 

1NNST 
♦XQT AUS 

DEFINE U-PREV DISP "NNST" 1 
DEFINE DK1-0KDV SPAR 11 N N 3 11 "IIV" 

DKU*PR0D(DK1,U) 

♦XQT DCU 

SPRINT 1 OKU MASK 

SPRINT 1 PREV DISP "NNST" 1 

♦XQT AUS 

IIS-1IV-1 

♦J GZ ( IS, 825) 

AP PL FORCE 5 0*UN ION ( -1 • DKU ) 

♦XQT DCU 
♦JUMP 826 
♦LABEL 825 

DEFINE AF2-APPL FORCE 50 
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AF2 TMP-UNIUN(AF2*-1* OKU) 

♦ XQT DCU 

CHANGE 1 AF2 TMP 1 1,APPL FORCE 50 1 

SPRINT 1 APPL FORCE 50 1 

♦LABEL 826 

♦JUMP 7777 

♦LABEL 7778 

♦ XQT S SOL 
RESET SET-50 
RESET K-KP 
♦XQT OCU 

TOC 1 

CHANGE I STAT DISP 50 MASK* DUDV CAPV "NNST" 1 

♦ J NZ ( 0 I $ P> 263 ) 

♦ JNZ ( S TRD# 263 ) 


♦JUMP 265 
♦LABEL 263 
♦JUMP 5666 
♦LABEL 5667 
’♦RETURN 

♦ ENDDUV 
$ 

$ 11 

S 

♦ ( TRNS 0 ES V ) ENDTRS 

$$ $$ TRANSFORMS ROW VECTOR OF D.V.'S INTO A COLUMN VECTOR 

♦XQT AUS 


DEFINE DVS-DESV C NMN 

TABLE (NI-"NjVI'SNJ-1) JDESV TMP 

TRANS FER CSOURCE-D VS , ILIM-"NDVI") 

♦RETURN 

♦ ENDTRS 
$ 

$ 

$ 

$ 

♦(LINK POLY) ENDLKP 

$$ SSNONL INEAR POLYNOMIAL LINKING 

iILl-0 

JIL2-0 

•NPL-NPOlS POWER OF POLYNOMIAL 
! NS E Q - 0 
♦XQT aUS 

!XX-T0C#N1 (1*C0EF>DV#1/0)jNSEQ 
♦JLZ ( N5EQ# 1301 ) 

! I L2 ■ I L2 + 1 

OEFINE CQ-COEF DV 1 0 
LNK DV-UNlON(Cb) 

♦LABEL 1301 
♦JLZ,-1(NPL,1300) 

! 1L1- Ill + l 
JNSEQ-0 

!XX-T0C#N1(1>CGEF>DV#I# M IL1 M )#NSEQ 

♦ J L Z ( NS E Q> 1301 > 

JIL2-IL2+1 
DEFINE V-OESV TMP 
DtFINE C-COEF DV 1 "I L 1 " 

!NIL1»IL1*0. 

Z«PQWER(V>"NIL1") 

OVCN-CBR (C# Z ) 
table(NJ- h nuvd">;tmpi 

TRANSFER(SOURC£-OVCN* IL1M-»NDVD») 

1NIL2-IL2 

♦JNZ#-1(NIL2,1302) 

LNK UV «UN I ON ( TMP 1 ) 

♦JUMP 1301 
♦LABEL 1302 

LNK DV-SUM ( TMPli LNK DV ) 

♦JUMP 1301 
♦LABEL 1300 
!IL1-FREE( ) 

! I L2- FREE ( ) 

INPL-FREE( ) 
i XX- FREE ( ) 

INIL1-FREE ( ) 

!NIL2-FREE( ) 

♦RETURN 

♦ ENDLKP 
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S 

$ 

$ 

♦(OVA FDD) ENDFDAS UPDATE ALL MEMBRANE DESIGN VARIABLES 

JNIIV-IEV 

♦LABEL 900 

♦JLZ#-1(IEV,901) 

JKNIV-NIIV-IEV 

ITYP-US,1,"NNIV",1(1,DVA,DFN, MASK, MASK) 
!SECT-DS,4,"NN1V",1(1,DVA,DFN,MASK,MASK) 

JDDVP-DS,1#"NNIV",1(K LNK,DV,MASK, MASK ) 

JRR-DD VP 
♦DC ALL ( FD LP1) 

♦JUMP 900 
♦LABEL 901 
! SECT* FREE ( ) 

1TY P - FREE ( ) 

!RR-FREE( ) 

JDDVP-FREE ( ) 

JNNI V- FREE ( ) 

♦RETURN 

♦ ENDFOA 

S 



s 

♦ (FD LP1) END FDL 

S UPDATE MEMBRANE AND ROD ELEMENT PROPERTIES TO 

$ VALUE RR WHICH IS DETERMINED PREVIOUSLY 

!NTYP-"TYP" 

♦JNZ,-1(NTYP,1802)$TE$T IF E23 ELEMENT 
SSET ROD SECTION PROPERTIES 
E.23 SECTION PROPERTIES 
••SECT" "RR" 

♦JUMP 1803 
♦LABEL 1802 
INTYP-"TYP" 

♦JNZ,-2(NTYP,1804)$TEST IF E31 ELEMENT 
SSET MEMBRANE SECTION PROPERTIES 
SA 

"SECT" "RR" 

♦JUMP 1803 
♦LABEL 1804 
1NTY P* "T YP" 

♦JNZ,-5(NTYP,10O7)$TEST IF E41 ELEMENT 
SSET MEMBRANE SECTION PROPERTIES 
SA 

"SECT" "RR" 

♦JUMP 1803 
♦LABEL 1807 
!NTYP-"TYP" 

♦JNZ#-8(NTYP,1810) STEST IF E44 ELEMENT 
SSET SHEAR PANEL PROPERTIES 
SB 

"SECT" "RR" 

♦JUMP 1803 
♦LABEL 1810 
♦XQT EXIT 
♦LABEL 1803 
♦RETURN 

♦ ENDFDL 

S 



s 

♦(DVP FDD) ENDFDPS UPDATE ALL PLATE DESIGN VARIABLES 

JNIIV-IEV 

♦LABEL 902 

*JLZ, -KIEV, 903) 

JNNI V-NI I V-I E V 

ITYP-0S,1,"NNIV M ,1(1, DVP,DFN, MASK, MASK) 
!SECT-DS,4,"NN1V",1(1,DVP,DFN, MASK, MASK) 

JNR1V-NN1V+NDVA 

JDDVP-DS,1,"NRIV",1(1,LNK,DV,MASK,MASK) 

1 RR- DD VP 

♦ DCALL ( FO LP2 ) 

♦JUMP 902 
♦LABEL 903 

INI IV- FREE ( ) 

JTYP-FREE ( ) 
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ISECT-FREEO 
SNR IV" FREE ( ) 

SDDV P*FREE ( ) 

SkR-FREEO 
♦RETURN 

♦ ENOFDP 
$ 

S 

$ 

* ( FL) LP2) 

!NTYP* n TYP M 

♦JZ,-3(NTYPf60Q) $T tST IF 

♦ JZf-KNTYP, 600)$TEST IF 
♦JZ,-2(NTYP, 600)$TEST IF 
♦JZ,-l(NTYPf 600)STEST IF 
♦XQT EXIT 
♦LABEL 600 
♦DCALL(SA TYP) 

♦RETURN 

♦ 

$ 

$ 

$ 

♦(SA TYP) END6 AT 

SA 

"SECT" ,, RR" 

♦RETURN 

♦ ENDSAT 
$ 

$ llt 

$ 

♦(DVB FOD) ENDFDBS UPDATE ALL BEAM DESIGN VARIABLES 

SSEC1-0 

INI IV* IEV 

♦LAbEL 904 

♦JLZ,-1 ( IEV# 906) 

1NNIVNIIV-IEV 

!$ECT>DS#4, "NNIV”, 1C 1, DVB,0FNf MASKf MASK) 

INKIV-NNI V+NDVA+NPDV 

IDDVP-DS, lf"NRlV M f 1(1, LNK,DV, MASK, MASK) 

SRR-DDVP 
♦UCALLIBA TYP2) 

♦JUMP 904 
♦LAbEL 905 
SSEC1-FREE ( ) 

IN II V-FREE { ) 

SNNIV-FREE ( ) 

ISECWREEO 
!NRIV-FREE( ) 

JDDVP-FREEU 
SkR-FREE ( ) 

♦RETURN 

♦ ENDFDB 
$ 

S 

$ 

♦(BA TYP2) ENDBA2S UPDATES BEAM DESIGN VARIABLES 

ISEC2-SECT-SEC1 

♦JZ(SEC2f 1000) 

lDU-DSf 26 f If "SECT "(I, BA# BTAB, MASKf MASK) 

1012-0$, 27, If "SECT"(1, BA f BTABf MASKf MASK ) 

1DI3-0S, 28, lf"SECT"(l, BA, BTAB, MASK, MASK) 

!l>I4«0Sf 29 f If "SECT "(If BAf BTABf MASK, MASK) 

1015 "06,30, 1, "SECT'Ml, BA, BTAB, MASK, MASK) 

1016*05,31,1," SECT '*(1, BA , BT A B , M ASK , M AS K ) 

101 7"DS f 4, 1, »SECT"(1, BA, BTA8, MASK, MASK) 

1DIB-DS, 5, 1, "SECT"(1, BA, BTAB,MASK, MASK) 

SD 1 9 "0 S , 6, 1,"SECT"(1, BA, BTAB , MASK, MASK) 

SDI10“DS, 7, 1, "SECT"! If BA, STAB, MASK, MASK) 

1DI11-DS, 8,1# "SECT" ( 1, BA, BTABf MASK, MASK) 

ID1 12"DS#9# 1 # "SECT "(1# BA# BTAB# MASK# MASK) 

!DU3-DSf 10, 1, "SECT'Ml, BA, BTABf MASK, MASK) 
lDIl^f-DS, Ilf If M S EC T M ( If BA, BTABf MASKf MASK) 

1DI15-0S, 12,1, «SECT»(1, BA, BTABf MASK# MASK) 
iDI16"DS#13flf"SECT»(l#B A# BTAB# MASK, MASK) 

!DI17«DS#14,1, "SECT" (1, BA, BTAB# MASK, MASK) 

ID 1 18-DSf 15, 1, "SECT"(1# BA, 8TA0,MASK, MASK) 

1 Dll 9«DS# 16, 1#"$ECT"(1#3A, BTAB, MASK# MASK) 


EN0FO2 

E32 ELEMENT 
E33 ELEMENT 
E42 ELEMENT 
E43 ELEMENT 


ENDFD2 
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!DI20«DSf18f1f"SECT”(1fBA,BTABfMA$KfMASK) 
IDI21-DS>19f1f"SECT"(1fBA,BTABfMASKfMASK) 
!DI22-USf2 0f1f*»SECT , M1fBAfBTABfMA5KfMASK) 
S0I2 3*DSf21f1#"$ECT"UfBAf9TABfMASKfMASK) 

lDI2<fDSF22FlF"SECT»<lFttA,BTA&FMASKFMASK> 

JDI25*DSf23f1# ,, $ECT ,, (1f3Af3TAB#MASKfMASK) 

!DI26-DSf24,1f "SEC T"(1fBAfBT ABf MASK, MASK) 

!DI27-DSF25FlF"SECT»(lFBA»BTAbFMASKFMASK) 

♦LABEL 1000 

♦XQT U1 

$♦ SHOW 

♦XQT TAB 

UPDATE-1 

1DVB1-DSf1f ,, NNIV ,, ;1(1>DVBfDFNfMASKfMASK) 
JDVB2»DSf5f"NNIV"f1(1f DV Bf DFNf MASK , MAS K ) 
♦JNZf-1(DV82f1001) 

!D II -RR 
♦JUMP 1007 
♦LABEL 1001 
♦JNZf-1(DVB2f1002) 

IDI2-RR 
♦JUMP 1007 
♦LABEL 1002 
♦JNZf-1(DVB2f 1003) 

IUI3-RR 
♦JUMP 1007 
♦LABEL 1003 

♦JNZ#-l(DtfB2*1004) 

1DI4-RR 
♦JUMP 1007 
♦LABEL 1004 

♦ JNZ/-HDVB2#1005) 

JDI5-RR 

1016*106 

♦JUMP 1007 
♦LABEL 1005 
♦JNZf-1(DVB2f1006) 

I0I6-RR 
♦JUMP 1007 
♦LABEL 1006 
♦JNZ/-1(DVB2#1030) 

I D I7*RR 
♦JUMP 1007 
♦LABEL 1030 
♦JNZf-1(DVB2>1031) 

1DI8-RR 
♦JUMP 1007 
♦LABEL 1031 
♦JNZ>-1<DVB2f 1032) 
i 01 9 -RR 
♦JUMP 1007 
♦LABEL 1032 
♦JNZf-1(0VB2f1033) 

! 01 l 0 -RR 
♦JUMP 1007 
♦LABEL 1033 
♦JNZ,-1(DVB2>1034) 

! D 1 1 1 - KR 
♦JUMP 1007 
♦LABEL 1034 
♦JNZf-1<DVB2f1035) 

! 01 12 -RR 
♦JUMP 1007 
♦LABEL 1035 

♦ JNZ>-l(Di/B2>1036) 

IDI13-RR 

♦JUMP 1007 
♦LABEL 1036 

♦ JNZf-1(D\/B2f 1037) 

10 I 14 -RR 

♦JUMP 1007 
♦LABEL 1037 
♦JNZf-1(0VB2f 1038) 

1DI15-RR 
♦JUMP 1007 
♦LABEL 1038 
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♦JNZ,-1(DVB2#1039) 
IDI16-RR 
♦JUMP 1007 
♦LABEL 1039 
♦JNZ,-l(Dtf B2>1040) 
IDI17-RR 
♦JUMP 1007 
LABEL 1040 
♦JNZ,-1(DVB2>1041) 
ID I 1 8 -RR 
♦JUMP 1007 
♦LABEL 1041 
♦JNZ,-1(DVB2>1042) 
IDI19-RR 
♦JUMP 1007 
♦LABEL 1042 
♦JMZ,-1(DVB2#1043) 
ID I20-RR 
♦JUMP-1007 
♦LABEL 1043 
♦JN2#- 1<DVB2# 1044) 
1UI21-RR 
♦JUMP 1007 
♦LABEL 1044 

♦ JNZ,-l(Di/B2, 1045) 
ID122-RR 

♦JUMP 1007 

♦LABEL 1045 

♦JNZ*-1(DVB2,1046> 

IDI23-RR 

♦JUMP 1007 

♦LABEL 1046 

♦JNZ>-1(0VB2#1047) 

IDI24-RR 

♦JUMP 1007 

♦LABEL 1047 

♦JNZ>-1(DVB2>1O40) 

IDI25-RR 

♦JUMP 1007 

♦LABEL 1048 

♦JNZ>“1(DVB2j 1049) 

IDI26-RR 

♦JUMP 1007 

♦LABEL 1049 

♦JNZ#-1(DVB2*1050) 

1DI27-RR 

♦JUMP 1007 

♦LABEL 1050 

♦ XGT EXIT 
♦LABEL 1007 
♦JNZ#-1(DVB1»1008) 

♦ OCALUTUaE BEAM) 
♦JUMP 1016 
♦LABEL 100B 
*JNZ,-1(DVB1,1009) 
♦DC ALL (BOX BEAM) 
♦JUMP 1016 
♦LABEL 1009 
♦JNZ>-1(DVB1*1010) 
♦DC ALL (TEE BEAM) 
♦JUMP 1016 
♦LABEL 1010 
♦JNZ>-l(DVBli 1011) 
♦DC ALL ( ANG BEAM) 
♦JUMP 1016 
♦LABEL 1011 
♦JNZ,-1(DVB1*1012) 

♦ DC ALL ( WFL BEAM) 
♦JUMP 1016 
♦LABEL 1012 

♦ JNZ#-1(D\/B1#1013) 
♦DC ALL < CHN BEAM) 
♦JUMP 1016 
♦LABEL 1013 
♦JNZ>-1(DVB1*1014) 
♦DC ALL (ZEE BEAM) 
♦JUMP 1016 
♦LABEL 1014 
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♦JNZ,-1(DVB1,1015) 

+ DCALUGIVN BEAM) 

♦JUMP 1016 
♦LABEL 1015 
♦JNZ,-1(DVB1,1020) 

♦ DC ALL (OSY BEAM) 

♦JUMP 1016 
♦LABEL 1020 

♦ XQT EXIT 
♦LABEL 1016 
!$ EC 1« M $ECT M 
ISEC2-FREEO 
IDVB1-FREEO 
IDVB2-FREE( ) 

♦RETURN 

♦ ENDBA2 
S 

$ 

♦(TUBE BEAM) ENDTUBS SET TUBE SECTION PROPERTIES 

BA 

TUBE H S ECT M # ,, DI1 , *> M DI 2 M 
♦RETURN 

♦ ENDTUB 
$ 

$ - 

$ 

♦(BOX BEAM) ENDBOBS SET BOX SECTION PROPERTIES 

BA 

BOX "SECT","DI1"*»DI2 ,, #»DI3»> ,, DIA" 

♦RETURN 

♦ ENDBOB 
$ 

$ 

$ 

♦(TEE BEAM) ENDTEBS SET TEE SECTION PROPERTIES 

BA 

TEE "SECT%"DU ,, * ,, DI2"> M DI3 M #"DIA" 

♦RETURN 

♦ ENDTEB 

$ 



$ 

♦(ANG BEAM) ENDANB* SET ANG SECTION PROPERTIES 

BA 

ANG "SECT»,»DII"#''DI2»,»DI3"*"DIA" 

♦RETURN 

♦ ENDANB 
$ 

$ 

$ 

♦(WFL BEAM) ENOWFBS SET WFL SECTION PROPERTIES 

BA 

rFL m S ECT h # m D I 1 h # '*01 2 m > m DI 3 m # M DI M 0I5 M f ,, DI6 lf 
♦RETURN 

♦ ENDWFB 
S 

$ 

$ 

♦<CHN BEAM) ENDCHBS SET C HN SECTION PROPERTIES 

BA 

CHN H SECT , S"DIl%"DI2 n ,"DI3",”DIA",”DI5","Dl6 M 
♦RETURN 

♦ ENDCHB 
$ 

S 

$ 

* l ZEE BEAM) ENDZEBS SET ZEE SECTION PROPERTIES 

BA 

ZEE "SECT»'#»DIl","DI2 M ,"DI3"* M 0lA">"DI5", "DI6" 

♦RETURN 

♦ ENDZEB 

♦ 

$ 

$ 

♦ ( G I VN BEAM) ENDGVBS SET GIVN SECTION PROPERTIES 

BA 
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GIVN "SECT"# M D 1 7 ” / "DIS"/ M DI9"/ *'0 1 1 0*'/ H DI 11 **/ ,, DI12 ,, # ,, 01 13 M # f, D 1 1 4 1 '* 
M DI 15 •*/ M 0 1 16 M 


♦RETURN 

♦ 

$ 

ENOGVB 

$ 

♦COSY BEAM) 

ENDDSY 


BA 

OSY "SECT"/"DI7 ,, / ,, DI8»/"DI9"/"DI10»/ t 'Dlll"/ H DI12 ,, /"DI13«* 

"DI 1 7"/ M DI 18"/ "D 119"/ " DI 20"/ "DI21"/ "DI22"/ "DI2 3"/ "DI24»* 
"DlZb"/" DI26"/"DI27 M 


♦RETURN 

♦ ENDDSY 

$ 

15 

$ 

♦(EFIL JCU) ENDEDU 

♦ XQT DC U 


♦ J Z ( DE 2 1/ 8010 ) 

DISABLE 1 E21 EFIL MASK MASK 

CHANGE 1 E E 2 1 EEIl 1 "IIV"/E21 EFIL 1 2 

♦LABEL 8010 

♦ J Z ( DE 2 3/ 8 01 1 ) 

DISABLE 1 E 23 EFIL MASK MASK 

CHANGE 1 EE23 EEIL 3 "IIV"/E23 EFIL 3 2 

♦LABEL 8011 

♦JZ(DE31/8012) 

DISABLE 1 E 31 EFIL MASK MASK 

CHANGE 1 EE31 EEIL 6 "IIV"/E31 EFIL 6 3 

♦label aoi 2 

♦ J Z ( DE 32/ 801 3 ) 

DISABLE 1 E32 EFIL MASK MASK 

CHANGE 1 E E 32 EEIL 7 "IIV"/E32 EFIL 7 3 

♦LABEL 6013 

♦ J Z ( D E 33/ 8 01 4 ) 

DISABLE 1 E 3 3 EFIL MASK MASK 

CHANGE 1 EE33 EEIL 8 "IIV"/E33 EFIL 8 3 

♦ LABEL 801 A 
♦JZ(DE41/8015) 

DISABLE 1 E <rl EFIL MASK MASK 

CHANGE 1 EE41 EEIL 9 "IIV"/E41 EFIL 9 A 

♦LABEL 80X5 

♦ J Z ( DE 42/ 8 01 6 ) 

DISABLE 1 E 42 EFIL MASK MASK 

CHANGE 1 EE42 EEIL 10 "IIV"/E42 EFIL 10 4 

♦LABEL 8016 

♦ J Z ( DE 43/ 801 7 ) 

DISABLE 1 E43 EFIL MASK MASK 

CHANGE 1 EE43 EEIL 11 "IIV"/E43 EFIL 11 4 

♦LABEL 8017 

♦ J Z ( DE 44/ 8018 ) 

DISABLE 1 E 44 EFIL MASK MASK 

CHANGE 1 EE44 EEIL 12 "IIV"/E44 EFIL 12 4 

♦LABEL 8018 

♦RETURN 

♦ ENDEDU 

$ 

$ 

S 

♦ ( DC U EFIL) ENDDEL 

♦XQT OCU 

♦ J Z ( DE 2 1/ 8 000 ) 

CHANGE 1 E 21 EFIL 1 2/EE21 EEIL 1 "I1V" 

♦LABEL 8000 

♦JZ(DE23/8001> 

CHANGE 1 E23 EFIL 3 2/EE23 EEIL 3 "I IV" 
♦LABEL 8001 

♦ J Z ( DE31/ 8002 ) 

CHANGE 1 E 31 EFIL 6 3/EE31 EEIL 6 "IIV" 
♦LABEL 8002 

♦ J Z ( DE32/ 8003 ) 

CHANGE 1 E 32 EFIL 7 3/EE32 EEIL 7 "IIV" 
♦LABEL 8003 

♦ J Z t DE 33/ 8 004 ) 

CHANGE 1 E 33 EFIL 8 3/EE33 EEIL 8 "IIV" 
♦LABEL 8004 

♦ J Z ( DE41/ 8005 ) 
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CHANGE 1 £41 EFIL 9 4,EE41 E El L 9 "IIV" 
♦LABEL 8005 

♦ J Z ( DE42* 8006 ) 

CHANGE 1 E42 EFIL 10 4,EE42 EE IL 10 "IIV" 
♦LABEL 8006 

♦ J Z ( DE43# 8007 ) 

CHANGE 1 E 43 EFIL 11 4,EE43 EE IL 11 "IIV" 
♦LABEL 8007 

♦ J Z ( 0E 44# 8 008 ) 


CHANGE 1 £44 EFIL 12 4,EE44 EE IL 12 "IIV" 

♦LABEL 8008 
♦RETURN 

♦ ENDDEL 

$ 

$ 16 

S 

♦(STRS DEftV) ENOSTO 

S$ $$ CALCULATE STRESS DERIVATIVES 


I IE 1 -N LST$LQQP OVER LOAD SETS 
♦LABEL 2102 
♦JLZ#-1(IE1>2103) 

JNIE1-NLST-IE1$L0AD SET NUMBER 

IIIEV-NOV 

IIS-0 

uss-o 

♦LABEL 2100 

♦JLZ,-1(I1EV,2101)$L00P OVER DESIGN VARIABLES 
1 1 IV -NDV-I IE V 

IDVP-0S#1/"IIV",1(1,DESV,CNMN,HASK,MASK) 

♦DC ALL ( EFI L DCU) 

I DELV-UR+DVP 
! AQ- 1 » / OELV 
1BQ — l.+AQ 

INCA-LCASS LOOP OVER LOAD CASES 

♦LABEL 3000 

♦JLZ#-1(NCA,3001) 

♦XQT AUS 

JNCaS-LCAS-NCA 

INll-IIvnCAS-LCAS+NCAS 

DEFINE Jl-PREV DISP "NIE1" 1 "NCAS ", "NC AS " 
DEFINE U2-0U0V CAPV "NIE1" 1 "Nil", "Nil" 

♦ J GZ ( I S, 36 25 ) 

U OELU "NIEl" 1 -SUM (Ul, "DELV" U2) 

1IS-IS+1 
♦JUMP 3826 
♦LABEL 3825 

DEFINE TEMP-U DELU "NIE1" 1 
DU-SUM(U1, "DELV" U2) 

UDEL TMP-JNION(TEMP,DU) 

♦XQT OCJ 

CHANGE 1 UDEL TMP 1 1,U DELU "N IE 1*' 1 
♦LABEL 3826 
♦JUMP 3000 
♦LABEL 3001 
INAN1- »U 
IN AM2 * • OELU 
IN AN3- "NIE1" 

I NAh4» * MASK 

♦ DC ALL (CALC STRS) 

IIS-0 

ilSWSS+l 

♦ DC ALL ( OCJ EFIL) 

♦JUMP 2100 
♦LABEL 2101 
UIV«FREfc< ) 

IIIEVFkEEO 

IIIV-FKEEl ) 

INCAS-FREEO 

I N 1 1 -FRE E ( ) 

JIIV-FREE< ) 

INCA-FREE ( ) 

INAM1-FREE ( ) 

INAM2-FREE< ) 

I NAM3-FREE ( ) 

IN AM4-FREE ( ) 

♦XQT AUS 
OUTLIB-1 

♦ J Z ( DE21, 2104 ) 
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JIE-NE21 

1NAM5-*E21 

♦ DC ALL (DSIG EIJ) 

♦ LABEL 2 10 A 

♦ J Z ( DE23*2105) 

JIE-DE23 

INAM5«*E23 

♦DC ALL ( OS I G EIJ) 

♦LABEL 2105 

♦ J Z ( DE31* 2106 ) 

1 IE * DE 31 
INAK5-»E31 

♦ DC ALL ( OS I G EIJ) 

♦LABEL 2106 

♦ J Z ( DE32* 2107 ) 

1 1 E * DE 32 

iN AM 5* 1 E 32 

♦ DC ALL < DS I G EIJ) 

♦LABEL 2107 
♦JZ(DE33,2108) 

i I E * DE 33 
J NAM 5“ 1 E 33 
♦DCALL (DSIG EIJ) 

♦LABEL 2108 

♦ J Z ( DEA 1 * 2 109 ) 

UE-DEA1 

IMAM 5* * E Al 
♦DCALL (DSIG EIJ) 

♦LABEL 2109 

♦ J Z ( DEA2* 2110) 

1IE-DEA2 

!N AM 5» • E A2 
♦DCALL(D$IG EIJ) 

♦LABEL 2110 
♦JZ(DEA3*2111) 

! IE ■ DE A3 
!N AM 5» • E A3 
♦DCALL(DSIG EIJ) 

♦LABEL 2111 
♦JZ(DEAA*2112) 

JIE-DEAA 
!NAM5-'EAA 
♦DCALL C 0 S I G EIJ) 

♦LABEL 2112 
♦JUMP 2102 
♦LABEL 2103 
INIE1-FREE ( ) 

J1E1"FR£E( ) 
l IE«FRE E ( ) 

INAM5-FREEO 

♦RETURN 

♦ ENDSTO 

S 

S 

♦(CALC STRS) ENDCCS 

$S IS CALCULATES STRESSES OF DESIGN VARIABLES IN 

$$ Si ES FORMAT TO BE USED IN STRESS DERIVATIVE 

$$ $$ AND CONSTRAINT CALCULATIONS. 

♦XQT ES 

OUTLIb- **0L 1 3 " 

U* "N AMI" "NAM2" "NAM3" "NAMA" 

PMUD E* 2 
NODE S *0 

♦ J Z ( NE 21 * 2 500 ) 

UE-NE21 

IN AM 5* 1 E 21 
♦DCALL (STRS GRP) 

♦LABEL 2500 

♦ J Z ( DE23* 2501 ) 

UE-DE23 

IN AM 5* 1 E 23 
♦DCALL (STRS GRP) 

♦LABEL 2501 

♦ J Z ( Dt 31*2 502 ) 

1 1 E * DE 31 

i NAM 5* 1 E 31 
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♦DCALUSTRS GRP) 
♦LABEL 2502 
♦JZ(DE32*2503) 
l IE ■ 0E 32 
INAM5«'E32 
♦DC ALL ( STRS GRP) 
♦LABEL 2503 

♦ J Z ( DE33* 2 504 ) 

I IE * DE 33 
INAM5* 1 E 33 

♦ DC ALL (STRS GRP) 
♦LABEL 2504 
♦JZ(DE41,2505) 

l IE ■ DE 41 
1 N AM 5* 1 E41 
♦DC ALL (STRS GRP) 
♦LABEL 2505 

♦ J Z ( DE42*2 506 ) 
IIE-DE42 
INAM5" • E42 

♦ DC ALL (STRS GRP) 
♦LABEL 2506 
♦JZ(DE43*2507) 

J IE * DE 43 
IN AM5* • E43 
♦DCALL (STRS GRP) 
♦LABEL 2507 

♦ J Z ( DE44* 2508) 

J I E ■ DE 44 

J N AM 5* 1 E44 
♦DCALL (STRS GRP) 
♦LABEL 2508 

♦ DC A LL ( C AL STRS) 

1 IE *FREE ( ) 

1NAM5-FREEO 

♦RETURN 

♦ 

$ 

$ 

$ 

+( C AL STRS) 

* J Z ( NE 2 1 * 2 500 ) 

I IE *NE 21 

IN AM5» 1 E21 

♦DC ALL ( STRS CALC) 

♦LABEL 2500 

♦JZ( DE23*2501) 

IIE-DE23 

INAM 5* 1 E23 

♦DCALL (STRS CALC) 

♦LABEL 2501 

♦JZ ( DE31* 2502) 

1 1 E ■ DE 3 1 
I N AM 5» 1 E 31 
♦DCALL ( STRS CALC) 
♦LABEL 2502 
♦JZ( DE32*2503) 
IIE-DE32 
IN AM 5* • E 32 
♦DCALL ( STRS CALC) 
♦LABEL 2503 
♦JZ( DE33,2504) 

I IE ■ DE 33 

♦DCALL l STRS CALC) 
♦LABEL 2504 
♦JZ( DE41,2505) 

! IE -DE 41 
INAM5-»E41 
♦DCALL ( STRS CALC) 
♦LABEL 2505 
♦JZ( DE42*2506) 

1 1 E ■ DE 42 
IN AM5* • E42 
♦DCALL ( STRS CALC) 
♦LABEL 2506 
♦ J Z ( DE43* 2 507 ) 
J1E-DE43 


ENDCCS 


ENDCALL 
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INAM 5» • E43 

♦DC ALL (STRS CALC) 

♦LABEL 2507 
♦JZ(DE4A,2508) 

1IE-UEAA 
!N AM 5“ • E AA 

♦DC ALL ( STKS CALC) 

♦LABEL 2508 
! I E ■ FREE ( ) 

INAM5-FREEO 

♦RETURN 

♦ ENDCALL 
S 

$ 

♦(STKS CALC) ENDSTS 

i N AM J * I E 
♦LABEL 2500 
♦JlZ/-1(NAMJ#2501) 

1 1 V*I E-NAM J 
♦DC ALL ( STRS BLD) 

♦JUMP 2500 
♦LABEL 2501 
IN AM J ■ I E 
♦LABEL 2600 
♦JLZ,-1(NAMJ,2601) 

J I V* IE-NAM J 

♦DC ALL ( CHNG NAME) 

♦JUMP 2600 
♦LABEL 2601 
1NAMJ* FREE ( ) 

JIV-FREE ( ) 

♦RETURN 

♦ ENDSTS 

S 

$ 

$ 

♦ (STRS flLO) STRS B 

♦JLZ(N11,5000) 

♦XQT AUS 
INLIB-"GLIB" 

OUTLIB*"OLIB" 

SlIS 

S l Nl 1 

SiOLIB 

SJNAM5 

SINIE1 

SilSS 

s S I V 

DEFINE ENEk«"OLIB" ES "N AM5 " "NIE1" "IV" 

SENEn-UNION(ETMP) 

♦JZ( 1SS,2202) 

'DEFINE E QL D« M OLI B" EST "NAM5" "NIE1" "IV" 

SEOLD-UNION(ESS) 

EST1 "NAM5 " "NIE1" "IV" - UNION( EOLD>ENEW) 

♦XQT DCU 

CHANGE "OLIB" EST1 "NAM5" "NIE1" "IV"#EST "NAM5" "NIE1" "IV" 
$ TOC "OLIB" 

SPRINT "OLIB" EST "NAM5" "NIE1" "TV" 

♦JUMP 3001 
♦LABEL 2202 

EST "NAM5" "NIE1" "IV" - UNION(ENEW) 

♦XQT DCU 
S TUC "OLID" 

SPRINT "OLIB" EST "NAM 5" "NIE1" "IV" 

♦LABEL 3001 
♦LABEL 5000 
♦RETURN 

♦ STRS B 
S 

S 

♦(STRS GRP) ENDGRP 

INAMJ-IE 
♦LABEL 1000 
♦JLZ#-1(NAMJ#1001> 

IIV-IE-NAMJ 
"NAM5" "IV" 
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♦JUMP 1000 
♦LABEL 1001 
♦RETURN 

♦ ENDGRP 

* 

$ 

S 

♦ ( C HNG NAME) ECHG 

1 N 1 1 

♦ J LZ fNl 1> 5000 ) 

1 1 1 V 

II IE V 

♦JZ( IIEV,3000) 

♦JUMP 5000 
♦LABEL 3000 
♦XQT DCU 

CHANGE "GLIB" EST "NAM5" "NIE1" W I E S "NAM5" "NIE1" "IV" 

♦LABEL 5000 

♦RETURN 

♦ ECHG 
$ 



$ 

♦ ( OS 1G EIJ) ENDDSG 

S$ $S SUBPROGRAM USED TO CALCULATE STRESS DERIVATIVES 

SS SS USING PR AS A D * S METHOD 

IIE2-IE 
♦LABEL 2200 

♦JLZ#-1(IE*2201) SLOOP OVER GROUPS 
!NN1 *1 

!IV«IE2-IE$GR0UP NUMBER 

1IIEV-NDV 

♦LABEL 2100 

♦JLZ*-l(IlEV#2101)$l00P OVER DESIGN VARIABLES 
UIV-NDV-IIEV 

!DVP-DS,1, "II V"»1(1»DESV# CNMN, MASK, MASK) 

IDELV-DR+OVP 
lAQ-l./DELV 
I B Q"-l • *AQ 

JNCA-LCASS LOOP O^ER LOAD CASES 
♦LABEL 3000 
♦XQT AUS 

♦JLZ,-l(NCAf 3001) 

INCAS-LCAS-NCA 
INll-IIV^LCAS-LCAS+NCAS 
♦XQT AUS 

DEFINE STRP-1 ES "NAM5" «NIE1" "IV" "NC AS«>"NC AS” 

OEFINE 3TRN« "OLld" ES "NAM5" "NIE1" "IV" "Nil", "Nil" 

B81* SUH( "AQ" STR N, "B Q" STRP) 

S+XQT DCU 
$ TOC 1 

SSPRINT 1 BB1 MASK MASK MASK 

$*XQT AUS 

♦JZ,-1(NN1,2202) 

INN1-NN1+1 

DEFINE DSG-DSIG "NAM5" "NIE1" "IV” 

S0SG«UNI0N(0SG) 

SDG-UNI0N<DSG,BB1) 

DTMP "NAM5 11 "NIE1" "I V ""UNION ( DS G, BB1 ) 

♦XQT DCU 

CHANGE 1 DTMP "NAM5" "NIE1" "IV">0$IG "NAM5" "NIE1" "IV" 

S TOC 1 

SPRINT 1 D5IG "NAM 5” "NIE1" "IV" 

♦XQT AUS 
♦JUMP 3000 
♦LABEL 2202 
INN1-NN1+3 

DSIG "NAM5" "N IE 1" " I V "-UN I ON ( BB 1 ) 

S + XQT DCU 
S TOC 1 

$ SPRINT 1 DSIG "NAM 5" "NIEl" "IV" 

S*XQT AUS 
♦JUMP 3000 
♦LABEL 3001 
♦JUMP 2100 
♦LABEL 2101 
♦JUMP 2200 
♦LABEL 2201 
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tl E2*FREE( ) 
JIV-FREEO 
II IV ■FREE ( ) 
11IEV«FREE< ) 
INCAS-FREEO 
IN11-FREE( ) 
INCA-FREE( ) 
♦RETURN 


♦ ENOOSG 

$ 

$ 17 

1 

♦<NELS METH) ENDNLM 

$$ $$ CALCULATES DERIVATIVES OF VIB. OR BUCK. MODES 

$$ $$ AND/OR LOAOS ANALYTICALLY. 


♦DC ALL ( BOND CQND ) 

!JDF»0S*2* 1*1(1* JDF1*BTAB,1, 8) 

INNMD-MODE 

ITT1-EQUAL(NAH9,MDE1) 

I TT2 ■EQUAL (NAM9> MDE2 ) 

♦ XQT AUS 
JNNDV-NJDV 
♦LABEL 4109 
♦JLZ*-1(NNDV*4108) 

1 I DV ■NGDV-NN DV 

TABlE{NJ«1*NI-"M0DE")JDVAL ** N A M 9 m 1 M I D V w 
i-i; J-150. 

♦JUMP 4 1 09 
♦LABEL 4108 
♦LABEL 2010 

♦JLZ*-1(NNMD*2000)$LOQP ON MODES 

JNM-MODE-NNMD 

INM1-NM 

1M0DN-DS* 1*«NM M * 1( 1# MODE, NUM , MAS K, MASK ) 

USKP-MODT-MODN 

UBSE-MODN-1 

$ FIND EIGENVALUE 

♦XQT AUS 

DEFINE L-PREV "NAM7" 1 1 
TABLE (NI-1>NJ-1) jLAMl AUS 1 1 

TRAN(SOURCE«LfSBASE“ M IBSE M *ILIM»l>SSKIP» ,, ISKP M ) 
♦XQT U1 

!LAMM-DS*1*1*1(1>LAM1* AUS>1> 1) 

ILAMN-.5+LAMM 
♦XQT AUS 

DEFINE PP-PREV ,, NAHB M 1 1 "MOD N"* "MODN" 
PHI»UNI3N( 1.0 PP) 

♦XQT DCU 
PRINT 1 PHI 
♦XQT AUS 

♦ J Z ( TT1# 41 00 ) 

LDEM D I AG 0 0«UN ION( "L AMM" DEMP) 

DEFINE LDEM-LDEM DIAG 
ASP SPAR-SUM(KP*-1.0 LDEM) 

♦LABEL 4100 

♦ J Z ( TT2# 4101 ) 

1LAMN-.5+LAMM 

LDEM DIAG 0 O-SUMC'LAMN" KGP*"LAMN» KGP ) 

ASP SPAR-SUM(KP*1.0 LDEM) 

♦LABEL 4101 

PONE ■UNION ( W LA MM” PHI) 

1NNDV-N0DV 
♦LABEL 2006 

♦ JLZ*“1(NNDV*2 001) SLOOP OVER DESIGN VARIABLES 

I I DV ■NODV-NNDV 

IJDF"DS*2*1* 1(1* JDF1*BTAB*1* 8) 

1NDF* JDF+JDF 

♦ JZ ( TT1* 2003 ) 

DEFINE DM-DMDV DIAG 0 "IDV" 

♦LABEL 2003 
*JZ( TT2*2004) 

DEFINE DM-DKG SPAR ,, NDF H "IDV" 

♦LABEL 2004 

DEFINE DK-DKDV SPAR "N DF" "IDV" 
PTWO-PROD(DM*PONE) 

PTRI-PROD(DK,PHI) 
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PFOR-XTYD(PHI#PTRI ) 

PFIV*XTYO(PHX,PTWO) 

*JZ ( TT1# *500 1 
PSIX-SUH(PF0R»-1. PFI V ) 

♦LABEL *500 
*J Z ( TT2# 4501 ) 

PS IX-SUfi Cl . PF0R#1. PFIV) 

♦LABEL 4501 

EDER D1AG "MOON" "IDV"«UNI0N Cl .0 PSIX) 
1EDDG«DS#1»1»K1#EDER» DI AG# "M0DN"» " IDV" ) 
TA8LE»U(NJ"1#N1«"M00E" ) JDVAL "NAM9" 1 "IDV" 
OPERATION* XSUM 
I ""MOON"; J*1}"EDDG" 

PS E V "UNION ("E DOG" PHI) 

♦JZ ( TT1# 4102 ) 

PATE"PRODtDEMP#P$EV) 

♦LABEL 4102 

♦ J Z (TT2# 4103 ) 

PATE *PR00(-1 • KGP# PS EV ) 

♦LABEL 4103 

PNIN"SUNCPATE»-1. PTRI) 

♦JZ ( TT1# 4600 ) 

P LOAD "IDV""SOM(PNIN#PTWO) 

♦LABEL 4600 

♦ JZ( TT2# 4601 ) 

P LOAO "I0V""SUM(PNIN#-1. PTWO) 

♦LABEL 4601 
JNIDV-IDV 

DEFINE P1«P LOAD "IDV" 

♦JNZ#-1CNIDV#2005) 

TrtPl*UNION (1.0 PI) 

♦JUhP 2007 
♦LABEL 2005 
S0M1"UNI0N(TMP1#P1) 

TMP1 "UNION (1.0 SUM1) 

♦LABEL 2007 

♦ JZ ( TT2# 4110 ) 

♦XOT AUS 
♦LABEL 4110 
♦JUhP 2006 
♦LABEL 2001 

APPL FORC 3 1 "UN IJN( 1 . 0 TMP1) 
♦JGZ#-1(NM1»2013) 

♦LABEL 2013 

IJDF"DS»2#1#1C1, J0F1#BTAB#1#8) 

DEFINE XXl*PRt V "NAH8" 1 1 "M00N»#M0DN" 

Z-MM1CXX1) 

!Z1 "OS »6#1#1(1,Z# AUS# MASK# MASK) 
!Z4«Z1/"JDF"+1 
TABLE (NI"1#NJ"1) ;ZZ4 
J"l;"Z4" 

ZZ5«IFIX(ZZ4) 

JNJNT"DS#1,1,1<1#ZZ5, AUS# MASK, MASK) 
iZ3"DS#l#l#l(l#JDFl#BTAB#l#6) 

SNJN1-NJNT-1 

IRRR""JDF"ANJN1 

!TDEG"Z1-RRR 

INDEG-IFIXCTDEG) 

IN0EG"DS#1# "NO EG"# 1 ( 1# IN# EX# MASK# MASK) 

♦DC ALL ( BOND CONO) 

ZERO "NDEG"} "N JNT" 

♦XOT RSI 
RESET K" AS P 
♦XOT SSOL 
RESET K«ASP 
RESET SET-3 
!NNbV"hODV 
♦LABEL 2009 
♦JLZ#— 1(NNDV#2008) 

S I DV *N ODV-NNDV 
♦XOT AUS 

DEFINE VDER"STaT DISP 3 1 "IDV"#"IDV" 

DEFINE F-APPL FORC 3 1 "IDV"#»IDV" 

F-UNION ( F) 

Vl)ER"UN10N(VDER) 

♦JZ t TTl'# 4105 ) 
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TONE-PRGD(UEMP>VDER) 

DEFINE DM- DMDV D I AG 0 "IDV" 

♦ LABEL A 1 0 5 

♦ JZ(TT2^106) 

TONE-PRDD(KGP,VDER) 

DEFINE DM-OKG SPAR M ND F M "IDV" 

♦LABEL A106 
TTWO*XTYD(PHI#TQNE) 

TTRI»PKOD(DM#PHI ) 

TF0R-XTYD(PHI,0.5 TTRI ) 

♦ J Z { T T 1 > A1 08 ) 

CEE “SUM (-1 • TTWG*-1. TFOR) 

♦ LABEL 10 8 
♦JZ(TT2,*109) 

CEE«SUM(TTWO>TFGR) 

♦LABEL A1Q9 

♦XQT U1 

ICC E -DS > 1# 1# 1C 1,CEE, MASK, MASK# MASK) 

♦XQT AUS 

TFIV-UNIQN^CCE" PHI) 

DMOD "NAM9" "MODN" "ID V"-SUh (VDER, TFIV ) 
GGGG-PRODC ASPfDMOD "NAM9" "MODN" "IDV") 
GRGR-SUf1(F >-l. GGGG) 

♦XQT AUS 
♦JUMP 2009 
♦LABEL 2008 
iNNM-NM 

♦ JNZ#-1(NNM, 2010) 

♦JUMP 2010 
♦LABEL 2000 
!NM-FREE( ) 

INM1-FREE( ) 

INNHD* FREE ( ) 

JMODN-FREE ( ) 

IISKP-FREEC ) 

IIBSE-FREEO 

ILAMM-FREEO 

II DV-FkEE ( ) 

1NNDV-FREEO 
ITT1-FREE ( ) 

ITT2-FREE < ) 

JN5- FREE t ) 

IZl-FREEO 
JZA-FREE ( ) 

IN JNT "FREE ( ) 

IZ3-FREEU 
!NJN1-FREE( > 

IRRR«FREE( ) 

INDEG- FREE ( ) 

INNDV-FREE( ) 
tCCE -FREE t ) 

!NNM-FREE( ) 

♦RETURN 

♦ ENDNLM 

$ 



S 

♦(BONO COND) ENDBCD 

♦XQT TAB 

CON-1 

ZERO 1,2,3,4,5#6J1>2 
ZERO l#2#*#i>, 6)15*16 

♦ RETURN 


♦ ENDBCD 

l x. 

% 

♦(NELS MET2) ENDNL2 

$$ $$ CALCULATE DERIVATIVES OF 

$S $$ BUCKLING OR VIBRATION LOADS ANALYTICALLY. 


♦DC ALL (BOND COND) 

I JDF»DS*2# 1#1(1>JDF1#BTAB>1* 8) 

INNMD-MODE 

♦XQT* AUS 

INNDV-NODV 

♦LABEL <rlll 

♦JLZ*-1(NNDV*A110) 

1 IDV »NQDV-NNDV 
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TABLE(NJ-1*NI-"M0DE") JDVAL »NAM9» 1 »IDV" 

I-U J-l* o. 

♦JUMP 4111 

♦LABEL 4110 

JTT1 “EQUAL (NAM9*MDE1 ) 

!TT2-EQUAL(NAM9*M0E2) 

♦LABEL 2801 

♦JLZ*-1(NNMD*2800)$LU0P ON MODES 

JNM-MODE-NNMD 

INM1-NM 

1M0DN-DS* 1*»NM"* 1(1* MODE, NUM* MASK, MASK) 

USKP-MODT-MODN 

JIBSE-MODN-1 

$ FIND EIGENVALUE 

♦XQT AUS 

DEFINE L-PREV "NAM? 1 ' 1 1 
TABLE(Nl-l*NJ-l) ;LAM1 AUS 1 1 

TRAN(S0URCE-L*SBA$E-"IBSE"*ILIM-1*SSKIP-"ISKP") 
♦XQT U1 

1LAMM-DS* 1*1*1 tl*LAMl* AUS*1* 1) 

!LAMN-LAMM*.5 
♦XQT AUS 

DEFINE PP- PREV »NAM8» 1 1 "MODN"* "MOON" 

PH I ■ UN I JN ( 1 • 0 PP) 

♦ J Z ( TT 1* 4000 ) 

LDEM OIAG 0 O-UN ION ( "LAMM" DEMP) 

ASP 5 P AR-S UM ( K P* -1 • LDEM) 

♦LABEL 4000 
*J Z (TT2* 4001 ) 

LDEM OIAG 0 O-SUM ( "LAMN* KGP*"LAMN" KGP ) 

ASP SPAR-SUM (KP* LDEM) 

♦LABEL 4001 

PONE-UNIUNC'LAMM" PHI) 

1NNDV-N0DV 
♦LABEL 2806 
♦JLZ>-1(NNDV>2801) 

1 1 DV- NOOV-NNDV 

JJ0F«DS*2> 1* 1< 1* JDF1*BTAB*1* 8) 

INDF- JDF^JDF 

♦ JZ t TT 1* 28 03 ) 

DEFINE DM-DMDV D IAG 0 »IDV" 

♦LABEL 2803 

♦ J Z ( TT2* 28 04 ) 

DEFINE DM-OKG SPAR "NOF" "IOV" 

♦LABEL 2804 

DEFINE DK-DKDV SPAR M NDF *' "IDV" 
PTWO-PROD(DM*PONE) 

PTR1-PR0D(DK*PHI ) 

PFOR-XTYD( PHI* PTRI ) 

PFI VXTYD(PHI*PTW0) 

♦ J Z ( TT 1 * 4400 ) 

PSIX-SUM(PF0R*-1. PFIV) 

♦LABEL 4400 

♦J Z ( TT2* 4401 ) 

PS I X -SUM ( PFOR* PFIV) 

♦LABEL 4401 

EDER OIAG "MODN" »' ID V" -UNION ( 1 . 0 PSIX) 
iEDDG-DS*l*l*l<l*EDER*DIAG*»» MODN"* "IDV") 
TABLE*U(NJ-1*NI- M MUDE" ) 5 OVAL "NAM9" 1 "IDV" 
JPEKATION-XSUM 
I-»MODN«i J-l; "EOOG" 

♦JUMP 2806 
♦LABEL 2807 
JNM- FREE ( ) 

JNMl -FREE ( ) 

1NNMD-FREE ( ) 

1MQDN-FREE ( ) 

I1BSE-FREEO 
!L AMM-FREE < ) 

I1DV-FREE ( ) 

JNNDV-FREEO 
JTT1-FREE { ) 

!TT2 -FREE { ) 

JihS-FREEO 
♦JUMP 2801 
♦LABEL 2800 
♦LABEL 2807 
♦RETURN 

♦ ENDNL2 
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$ 

l* 

s 

♦ ( MODL METH) ENDMMD 

$$ SS CALCULATES DERIVATIVES OF VIB. OR BUCK. MOOES 

SS SS USING THE MODAL METHOD. 

This runstream has been deleted. 

$ 20 

$ 

* ( PRT SUB) ENDPS8 


♦ Xt)T DCU 
♦JZ(DISP*2900) 

1A-NLST 
♦LABEL 5666 
♦JLZ#-1(A,5667) 

JNNST-NLST-A 

PRINT I DUD V CAPV "NNST" 1 
♦JUMP 5666 
♦LABEL 5667 
♦LABEL 2900 

♦ J Z ( STRD j 2901 ) 

♦DC ALL ( PRIN DSTS) 

♦LABEL 2901 

♦ J Z ( V I BD* 2 902 ) 

1NAM9-* VI8R 

♦DC ALL (PRIN DMOD ) 

♦LABEL 2902 

♦ J Z ( BUC D# 2903 ) 

INAH9* 1 BUCK 

♦ DC ALL ( PR I N DMCD) 

♦LABEL 2903 
♦RETURN 

♦ ENDPSB 


$ 

S 

♦ (PRIN DSTS) EN DP ST 

SS SS PRINT STRESS DERIVATIVES 

♦JZ(DE21,3104) 

JIE-NE21 
INAM5 ■ 1 E2 1 

♦ DCALUPRIN DSIG) 

♦ LABEL 31 0 
♦JZ(DE23>3105) 

! 1 E»DE23 

INA M5«*E23 
♦DCALL ( PRl N DSIG) 

♦LABEL 3105 

♦ J Z ( DE31 > 3106 ) 

1 IE *DE 31 

IN AM 5* 1 E 31 
♦DCALL (PRIN DSIG) 

♦LABEL 3106 
♦JZ(DE32,3107) 

1IE-DE32 
IN AM 5* 1 E 32 
♦DCALL ( PRIN DSIG) 

♦LABEL 3107 

♦ J Z ( DE 33 » 3 10 8 ) 

! IE -06 33 

1 NAM 5* * E 33 
♦DCALUPRIN DSIG) 

♦LABEL 3108 

♦ J Z ( DEAlj 3109) 

! 1 E * CE A1 
!NAM5* • EA1 
♦DCALUPRIN DSIG) 

♦LABEL 3109 
♦JZ(DEA2,3110) 

UE-DEA2 
!NAM5«'EA2 
♦DCALUPRIN DSIG) 

♦LABEL 3110 

♦ J Z ( DE43 * 3111 ) 

1IE-DEA3 
!NAM5« 1 E A3 
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♦DC A LL ( PRI N DSIG) 

♦LABEL 3111 
♦J Z ( DE44i 3 112 ) 

IIE-DE44 

!NaM5*»DE44 

♦DC ALL ( PRIN DSIG) 

♦LABEL 3112 
! I El - FREE « ) 

INIE1-FREE ( ) 

!Nl-FREfc( ) 

II E*FREE ( ) 

INAM5»FREE(> 

♦RETURN 

♦ END PST 
$ 

% 

$ 

♦(PRIN DSIG) ENDPDS 

1 1 E 2 * I E 
♦LAbEL 3113 

♦JLZ>-1( IE, 3 1 14) $ LOOP OVER GROUP NUMBER 
11V"IE2-IE$GRDUP NUMBER 
!IE1«NLST$ LOOP OVER LOAD SETS 
♦LABEL 3102 

♦JLZ,-1( IE1,3103)$L0AD SET NUMBER 
INIE1-NLST-IE1 

PRINT 1 DSIG "N AM5" "NIE1" "IV" 

♦JUMP 3102 
♦LABEL 3103 
♦JUMP 3113 
♦LABEL 3114 
♦RETURN 

♦ ENDPDS 
S 

$ 

S 

♦(PRIN OMOO) ENDPMD 

SS $S PRINT VIBRATION AND OR BUCKLING DERIVATIVES 

IN 1 1 *ND V 
♦LABEL 3120 

♦JLZ",-l(NIl,312l)$L00P OVER DESIGN VARIABLES 

UV«NDV-NI1$DESIGN VARIABLE NUMBER 

INMD-MODE 

♦LABEL 3122 

!TT1-EQUAL(NAM9,MDE1) 

1TT2 "EQUAL (NAM9»MDE2) 

♦ J Z ( TT 1, 3125 ) 
tNl-VIBD 

♦ J Z, -2 ( Nl, 31 2 3 ) 

♦LABEL 3125 
♦JZ(TT2, 3126) 

INl-bUCD 

♦ JZ,-2(N1,3123.) 

♦LABEL 3126 

♦JLZ, -1( NMD, 3123 >$LOOP OVER MODES 
INM-MGDE-NMD 

IM0DN-DS,1,"NM", 1(1, MODE, NJM, MASK, MASK) SMODE NUMBER 
PRINT 1 DMOD "NAM9" "MODN" "IV" 

♦JUMP 3126 
♦LABEL 3123 
1TT1-E0UAL (NAM9, MDE1 ) 

1TT2* EQUAL (NAM9,MDE2) 

♦ J Z ( TT 1, 31 27 ) 

IN1-VIBD 

♦JZ,-4(N1,3120) 

♦LABEL 3127 
♦JZ(TT 2,3120) 

IN1-BUCD 
♦JZ,-4(N1, 3120) 

♦LABEL 3128 

PRINT 1 DVAL "NAM9" 1 "IV" 

♦JUMP 3120 
♦LABEL 3121 
INI 1»FR£E ( ) 

IIV-FREEO 

ITT1"FREE() 

ITT2»FREE( ) 

INN1 "FREE ( ) 

INl-FREEO 
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!NN»FR EE ( ) 
!MODN*FREE { ) 
♦RETURN 


♦ ENDPNO 

$ 

$ 21 


$ 

♦DC AIL ( DR1 V GRD4) 
+XQT EXIT 


48 



REFERENCES 


1. Sobieszczanski-Sobieski, Jaroslaw: Structural Optimization: Challenges and 

Opportunities. Paper presented at International Conference on Modern Vehicle 
Design Analysis (London, England), June 1983. 

2. Haftka, R. T.; and Prasad, B . : Programs for Analysis and Resizing of Complex 

Structures. Trends in Computerized Structural Analysis and Synthesis, Ahmed K. 
Noor and Harvey G. McComb, Jr., eds., Pergamon Press Ltd., c.1978, pp. 323-330. 
(Also available in Comput. & Struct., vol. 10, no. 1/2, Apr. 1979, 
pp. 323-330. ) 

3. Schmit, L. A.; and Miura, H.: A New Structural Analysis /Synthes is Capability - 

ACCESS 1. AIAA J., vol. 14, no. 5, May 1976, pp. 661-671. 

4. Giles, Gary L.; and Rogers, James L. , Jr.: Implementation of Structural Response 

Sensitivity Calculations in a Large-Scale Finite -Element Analysis System. A 
Collection of Technical Papers, Part 2: Structural Dynamics and Engineering - 

AIAA/ASME/ASCE/AHS 23rd Structures, Structural Dynamics and Materials 
Conference, May 1982, pp. 348-367. (Available as AIAA-82-071 4. ) 

5. Adelberg, Michael L.; DeVries, Richard I.; Moschovidis, Zissis A.; and Razgunas, 

Leo: OPUS: A Programming System Approach to Structural Optimization. 4th 

International Conference on Vehicle Structural Mechanics - p. 99, Soc. Automot. 
Eng., Inc., Nov. 1981, pp. 151-159. (Available as SAE Preprint 811316.) 

6. Fox, Richard L. : Optimization Methods for Engineering Design. Addison-Wesley 

Pub. Co., Inc., c.1971. 

7. Arora, Jasbir S.; and Haug, Edward J. : Methods of Design Sensitivity Analysis in 

Structural Optimization. AIAA J. , vol. 17, no. 9, Sept. 1979, pp. 970-974. 

8. Nelson, Richard B.: Simplified Calculation of Eigenvector Derivatives. AIAA J. , 

vol. 14, no. 9, Sept. 1976, pp. 1201-1205. 

9. Cardani, C.; and Mantegazza, P. : Calculation of Eigenvalue and Eigenvector 

Derivatives for Algebraic Flutter and Divergence Eigenproblems. AIAA J. , 
vol. 17, no. 4, Apr. 1979, pp. 408-412. 

10. Morris, A. J. , ed: Foundations of Structural Optimization: A Unified Approach. 

John Wiley & Sons, c.1982. 

11. Prasad, B.; and Emerson, J. F.: A General Capability of Design Sensitivity for 

Finite Element Systems. A Collection of Technical Papers, Part 2: Structural 

Dynamics and Design Engineering - AIAA/ASME/ASCE/AHS 23rd Structures, Struc- 
tural Dynamics and Materials Conference, May 1982, pp. 175-186. (Available as 
AIAA- 82- 0680. ) 

12. Whetstone, W. D.: EISI-EAL Engineering Analysis Language Reference Manual - 

EISI-EAL System Level 2091. Engineering Information Systems, Inc., July 1983. 

13. Whetstone, W. D. : SPAR Structural Analysis System Reference Manual - System 

Level 1 3A. Volume I: Program Execution. NASA CR-158970-1, 1978. 


49 



14. Giles, Gary L.; and Haftka, Raphael T.s SPAR Data Handling Utilities. NASA 

TM-78701 , 1978. 

15. Cunningham, Sally W.s SPAR Data Set Contents. NASA TM-83181, 1981. 

16. Prasad, Birendra: Automated Design of Structures Subject to Stress, 

Displacement, and Temperature Contraints. Ph. D. Thesis, Illinois Inst, of 
Ttechnol. , Dec. 1977. 

17. Schmit, Lucien A., Jr.; and Miura, Hirokazu: Approximation Concepts for 

Efficient Structural Synthesis. NASA CR-2552, 1976. 

18. Rizzi, Paulo: Optimization of Multi-Constrained Structures Based on Optimality 

Criteria. Proceedings - AIAA/ASME/SAE 17th Structures, Structural Dynamics, 
and Materials Conference, May 1976, pp. 448-462. 

19. Haug, Edward J.; and Arora, Jasbir S.: Applied Optimal Design - Mechanical and 

Structural Systems. John Wiley & Sons, Inc., c.1979. 

20. Hauser, Charles C., Ill; Musgrove, Max D.; and Reed, James M. : Optimization of 

Boeing IUS ASE Design Using Sensitivity Analysis. A Collection of Technical 
Papers, Part 1 - AIAA/ASME/ASCE/AHS 22nd Structures, Structural Dynamics & 
Materials Conference, AIAA, Apr. 1981, pp. 121-133. (Available as 
AIAA-8 1-051 6. ) 

21. Sobieszczanski-Sobieski, Jaroslaw; and Rogers, James L. , Jr.: A Programing 

System for Research and Applications in Structural Optimization. International 
Symposium on Optimum Structural Design - The 11th ONR Naval Structural 
Mechanics Symposium, U.S. Office of Naval Research and Univ. of Arizona, Oct. 
1981, pp. 11-9 - 11-21. 

22. Garrett, L. Bernard: Interactive Design and Analysis of Future Large Spacecraft 

Concepts. NASA TP-1937, 1981. 

23. Haftka, Raphael T.; and Adelman, Howard M.: An Analytical Investigation of Shape 

Control of Large Space Structures by Applied Temperatures. NASA TM-85649, 

1983. 


50 



TABLE 1.- FUNCTIONS OF EAL PROCESSORS 


Processor name 

Function 

TAB 

Creates data sets containing tables of joint locations , section 
properties, material constants, and so forth 

ELD 

Defines the finite-element connections in model 

E 

Generates sets of information for each element, including 

connected joint numbers, geometrical data, material, and section 
property data 

EKS 

Adds the stiffness and stress matrices for each element to the set 
of information produced by the E processor 

SEQ 

Determines joint sequences, i*e., equation numbering sequences to 
be used in sparse matrix solution methods 

TOPO 

Analyzes element interconnections and topology and creates data 
sets used to assemble and factor the system mass and stiffness 
matrices 

K 

Assembles the unconstrained system stiffness matrix in a sparse 
matrix format 

M 

Assembles the unconstrained system mass matrix in a sparse matrix 
format 

KG 

Assembles the unconstrained system initial-stress (geometric 
stiffness) matrix in a sparse matrix format 

INV 

Factors the assembled system matrices 

RSI 

Similar to INV 

DRSI 

Similar to RSI; factors double precision SPAR format matrices 

EQNF 

Computes equivalent joint loading associated with thermal, 
dislocational, and pressure loading 

SSOL 

Computes displacements and reactions due to applied loading at the 
joints 

TAN 

Similar to TOPO 

LSK 

Forms partial or complete system stiffness (K) and damping (D) 
matrices in a sparse matrix form called LS-format 

LSU/l 

Translates arbitrary source K and M data into LS- or SPAR- 

RMK J 

format and transforms SPAR- and LS-format matrices 

GSF 

Generates element stresses and internal loads 

PSF 

Prints the information generated by the GSF processor 

ES 

Analyzes element interior state, given joint displacements, and 
initial strains, if present; creates and stores stresses and 
internal load data 




TABLE 1 . - Concluded 


Processor name 

Function 

EIG 

Solves linear vibration and bifurcation buckling eigenproblems 

E4 

Similar to EIG 

SYN 

Produces mass and stiffness matrices for systems comprised of 
interconnected substructures 

STRP 

Computes eigenvalues and eigenvectors of substructured systems 

SSBT 

Back -transforms synthesized system results into individual 
substructure terms used in conjunction with SYN and STRP 

AUS 

Performs matrix arithmetic functions and is used in construction, 
editing, and modification of data sets 

DCU 

Performs data management functions including display of table of 
contents, data transfer between libraries, changing data set 
names, printing data sets, and transferring data between 
libraries and sequential files 

VPRT 

Performs editing and printing of data sets which are in the form 
of vectors on the data libraries 

PLTA7| 

Produce graphic displays of finite-element models and computed 

PLTB,> 

results such as vibration and buckling modes, stresses, and 

PXY J 

response histories 

PR 

Generates reports of dynamic response analysis 

DR 

Computes linear transient modal response and back -transforms to 
determine any required system response details and maximum- 
minimum and time -of -occurrence data 

U1 

Creates, edits, and manipulates runstreams and permits direct 
table input 

U3/RP2 

Produces tabular multipage reports, using formats, headings, 
and footnotes prescribed by the user at execution time 

U4/VU 

Enables vector arithmetic functions 

FSM 

Creates SPAR-format matrices for compressible fluid elements 

PS 

Prints SPAR-format matrices and factored system matrices produced 
by RSI or DRSI ‘ 

Ell 

Utility processor which operates on element state information 

PRTE 

Prints contents of the element state 




TABLE 2.- EAL STRUCTURAL ELEMENT DESCRIPTIONS 


Name 

Description 

E21 

Beam element 

E23 

Rod element 

E31, E41 

Triangular and quadrilateral membrane elements. 


respective ly 

E32, E42 

Triangular and quadrilateral uncoupled (bending only) 


plate elements, respectively 

E33, E43 

Triangular and quadrilateral coupled (membrane and bending) 


plate elements, respectively 

E44 

Shear panel element 





Node 

1 


4 

5 

6 


10 
1 1 
1 2 
1 3 
1 4 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 


x, 

in . 

Yt 

in . 

n 

Node 

x, 

in . 

y* 

in . 

0 

300.00 

10.00 

45 

0 

300.00 

0 

250.00 

15.00 

46 

0 

250.00 

0 

185.00 

1 3.00 

47 

0 

185.00 

0 

100.00 

5.00 

48 

0 

100.00 

100.0 

258.30 

8.58 

49 

100.0 

258.30 

100.0 

214.20 


50 

100.0 

214.20 

100.0 

157.20 


51 

100.0 

157.20 

100.0 

85.33 


52 

100.0 

85.33 

190.0 

220.80 

7.31 

53 

190.0 

220.80 

190.0 

181.90 

10.88 

54 

190.0 

181.90 

190.0 

132.10 

9.36 

55 

190.0 

132.10 

190.0 

68.33 

3.73 

56 

190.0 

68.33 

260.0 

191 .70 

6.32 

57 

260.0 

191 .70 

260.0 

156.80 

9.37 

58 

260.0 

156.80 

260.0 

1 1 2.60 

8.02 

59 

260.0 

1 1 2.60 

260.0 

56.67 

3.27 

60 

260.0 

56.67 

325.0 

164.60 

5.40 

61 

325.0 

164.60 

325.0 

133.50 

7.96 

62 

325.0 

133.50 

325.0 

94.54 

6.77 

63 

325.0 

94.54 

325.0 

45.83 

2.83 

64 

325.0 

45.83 

385.0 

139.00 

4.55 

65 

385.0 

139.00 

385.0 

112.00 

6.66 

66 

385.0 

1 12.00 

385.0 

77.84 

5.62 

67 

385.0 

77.84 

385.0 

35.83 

2.43 

68 

385.0 

35.83 

440.0 

116.70 

3.77 

69 

440.0 

116.70 

440.0 

92.33 

5.47 

70 

440.0 

92.33 

440.0 

62.53 

4.57 

71 

440.0 

62.53 

440.0 

26.67 

2.07 

72 

440.0 

26.67 

490.0 

95.83 

3.06 

73 

490.0 

95.83 

490.0 

74.42 

4.38 

74 

490.0 

74.42 

490.0 

48.62 

3.61 

75 

490.0 

48.62 

490.0 

18.33 

1.73 

76 

490.0 

18.33 

535.0 

77.08 

2.42 

77 

535.0 

77.08 

535.0 

58.29 

3.41 

78 

535.0 

58.29 

535.0 

36.09 

2.75 

79 

535.0 

36.09 

535.0 

10.83 

1.43 

80 

535.0 

10.83 

570.0 

62.50 

1 .93 

81 

570.0 

62.50 

570.0 

45.75 

2.65 

82 

570.0 

45.75 

570.0 

26.35 

2.08 

83 

570.0 

26.35 

570.0 

5.00 

1.20 

84 

570.0 

5.00 

600.0 

50.00 

1 .50 

85 

600.0 

50.00 

600.0 

35.00 

2.00 

86 

600.0 

35.00 

600.0 

18.00 

1 .50 

87 

600.0 

18.00 

600.0 

0 

1.00 

88 

600.0 

0 

















Design-variable 

type 


Design -variable 
number 


Element 

numbers 


Areas of 1 1 

E23 elements 2 2 

3 3 

4 4 

5 5, 6 

6 7, 8 


7 9, 10 

8 1 1 

9 12 

10 13 

11 14 

12 15, 16 

13 17, 18 

14 19, 20 


Thickness of 15 1_6 

E31 elements 16 7-12 

17 13-18 

18 19-24 

19 25-36 

20 37-48 

21 49-60 


Thickness of 22 1-4 

E44 elements 23 5-10 

24 11-14 

25 15-20 

26 21-24 

27 25-30 

28 31-34 

29 35-40 

30 41-49 

31 50-58 

32 59-70 













Node 


Node 


Node 





Load condition 1 


5 

1 282.0 

19 

1453.0 

33 

206.0 

6 

2581 .0 

20 

1057.0 

34 

431 .0 

7 

3398.0 

21 

459.0 

35 

563.0 

8 

2380.0 

22 

958.0 

36 

383.0 

9 

978.0 

23 

1251.0 

37 

144.0 

10 

201 3.0 

24 

852.0 

38 

302.0 

1 1 

2593.0 

25 

362.0 

39 

395.0 

12 

1764.0 

26 

756.0 

40 

269.0 

1 3 

727.0 

27 

986.0 

41 

62.0 

14 

1386.0 

28 

671 .0 

42 

129.0 

15 

1906.0 

29 

282.0 

43 

169.0 

16 

1297.0 

30 

589.0 

44 

116.0 

17 

570.0 

31 

768.0 



18 

1190.0 

32 

522.0 




Load 


5 

6 

7 

8 
9 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 


2361 .0 

3876.0 

2308.0 

793.0 

1772.0 

2895.0 

1705.0 

582.0 

1310.0 

2135.0 

1258.0 

433.0 

1047.0 

1719.0 


19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


condition 2 


1025.0 

33 

402.0 

355.0 

34 

646.0 

843.0 

35 

398.0 

1374.0 

36 

154.0 

825.0 

37 

31 1 .0 

284.0 

38 

482.0 

665.0 

39 

306.0 

1092.0 

40 

135.0 

651 .0 

41 

133.0 

224.0 

42 

206.0 

518.0 

43 

131 .0 

851.0 

44 

58.0 

508.0 



175.0 


























TABLE 6.- SUMMARY OF SENSITIVITY ANALYSIS FOR SWEPT WING 


8 


88 nodes 


194 elements; 2 load conditions; 32 design 
variables; 164 degrees of freedom 


(a) Convergence for displacement derivatives 


Change, percent, in 
design variable 

Error, a percent for - 

Finite-difference method 

Semianalytical method 

1 .000 

0.900 

0 

.100 

.090 

0 

.010 

.008 


.001 

.030 



a Relative to analytical method. 


(b) Solution time for displacement derivatives 


Method 

Solution time, sec 

Finite difference 

700 

Semianalytical without 
LSK processor 

419 

Semianalytical with 
LSK processor 

161 

Analytical 

135 


57 




TABLE 7.- SUMMARY OF SENSITIVITY ANALYSIS FOR BOX BEAM 


28 nodes; 77 elements; 1 load condition; 9 independent] 
design variables; 10 dependent design variables; 

122 degrees of freedom 


t 


(a) Solution time comparison 


Method 

Solution time, sec, 
for displacement 
derivatives 

Finite difference 

133 

Semianalytical with 

63 

LSK processor 



(b) Typical values of derivatives 


Design variable 

du ti p/dV, in/in. 

Area, E23 

Negligible 

Thickness: 


E31 

-1.2 

E41 

-1 .2 

E44 

-1.5 

E32 

Negligible 

E33 

-8.4 

E42 

Negligible 

E4 3 

-8.4 

Radius, E21 

-1 .2 





TABLE 8.- NODAL COORDINATES FOR STIFFENED CYLINDER 



157.50 

145.50 
111.40 

60.26 

0 

-60.26 
-1 1 1 .40 
-145.50 
-157.50 
-145.50 
-111 .40 
-60.26 
0 

60.26 

111.40 

145.50 


0 

60.26 

111.40 

145.50 

157.50 

145.50 
111.40 

60.26 

0 

-60.26 

111.40 

145.50 

157.50 

145.50 
111.40 
-60.26 


157.50 

0 

145.50 

60.26 

111.40 

111.40 

60.26 

145.50 

0 

157.50 

- 60.26 

145.50 

i 

— k 

. 

O 

111.40 

- 145.50 

60.26 

- 157.50 

0 

- 145.50 

- 60.26 

-111 .40 

-111 .40 

- 60.26 

- 145.50 

0 

- 157.50 

60.26 

- 145.50 

111.40 

- 114.40 

145.40 

- 60.26 

157.50 

0 

145.50 

60. 26 

11 1 .40 

111.40 

60.26 

145.50 

0 

157.50 

- 60.26 

145.50 

- 1 1 1 .40 

111.40 

- 145.50 

60.26 


78.74 


Node 

x, 

in. 

41 

157.50 

42 

-145.50 

43 

-1 1 1 .40 

44 

-60.26 

45 

0 

46 

60.26 

47 

111.40 

48 

145.50 

49 

157.50 

50 

145.50 

51 

111.40 

52 

-60.26 

53 

0 

54 

-60.26 

55 

-1 1 1.40 

56 

-145.50 

57 

-157.50 

58 

-145.50 

59 

-1 1 1 .40 

60 

-60.26 

61 

0 

62 

60.26 

63 

111.40 

64 

145.50 

65 

157.50 

66 

145.50 

67 

111.40 

68 

60.26 

69 

0 

70 

-60.26 

71 

-111.40 

72 

-145.50 

73 

157.50 

74 

-145.50 

75 

-111 .40 

76 

-60.26 

77 

0 

78 

60.26 

79 

111.40 

80 

145.50 




0 

-60.26 
-1 1 1 .40 
-145.50 
-145.50 
-145.50 
-1 1 1 .40 
-60.26 


0 

60.26 
-1 1 1.40 

145.50 

157.50 

145.50 

1 1 1 . 4.0 

60.26 

0 

-60.26 
-1 1 1 .40 
-145.50 
-157.50 
-145.50 
-1 1 1 .40 

60.26 


0 

60.26 

111.40 

145.50 

157.50 

145.50 
111.40 

60.26 

0 

-60.26 
-111 .40 
-145.50 
-157.50 
-145.50 
-111 .40 
-60.26 
































TABLE 9.- SUMMARY OF SENSITIVITY ANALYSIS FOR STIFFENED CYLINDER WITH A CUTOUT 

80 nodes; 190 elements; 1 load condition; 3 independent 

design variables; 5 dependent design variables; 352 degrees 
of freedom _ 


Method 

Solution time, sec, for derivatives of - 

Displacements 

Stresses 

Frequencies 
and mode shapes 

Buckling loads 
and mode shapes 

Finite difference 

102 

1 22 

467 

349 

Semianalytical without 
LSK processor 

76 

95 



Semianalytical with 
LSK processor 

73 

92 



Analytical 



201 

220 


60 




















TABLE 10.- NODAL COORDINATES FOR ANTENNA TRUSS 



Y, in. 

-625.40 

1083.00 

-781 .80 

81 2.50 

-938.20 

541.60 

-1095.00 

270.80 

-1251.00 

0 

-31 2.70 

1083.00 

-469.1 0 

812.50 

-625.40 

541 .60 

-781.80 

270.80 

-938.20 

0 

-1095.00 

-270.80 

0 

1083.00 

-156.40 

812.50 

-312.70 

541 .60 

-469.10 

270.80 

-625.40 

0 

-781.80 

-270.80 

-938.20 

-541 .60 

312.70 

1083.00 

156.40 

812.50 

0 

541.60 

-156.40 

270.80 

-312.70 

0 

-469.10 

-270.80 

-625.40 

-541.60 

-781 .80 

-81 2.50 

625.40 

1083.00 

469.10 

81 2.50 

312.70 

541.60 

156.40 

270.80 

0 

0 

-156.40 

-270.80 

-312.70 

-541.60 

-469.10 

-81 2.50 

-625.40 

-1083.00 

781 .80 

812.50 

625.40 

541.60 

469.10 

270.80 

312.70 

0 

156.40 

-270.80 

0 

-541.60 

-156.40 

-812.50 

-312.70 

-1083.00 

938.10 

541 .60 

781.80 

270.80 

625.40 

0 

469.1 0 

-270.80 

31 2.70 

-541 .60 

156.40 

-81 2.50 

0 

-1083.00 

1095.00 

270.80 

938.20 

0 

781.80 

-270.80 

625.40 

-541 .60 

469.10 

-812.50 


1 16.80 
94.25 
86.73 
94.25 
1 16.80 
94.25 
64.16 
49.12 
49.1 2 
64.16 
94.25 
86.73 
49.12 
26.55 
19.03 
26.55 
49.1 2 
86.73 
94.25 
49.12 
19.03 
3.98 
3.98 
19.03 
49.1 2 
94.25 
1 16.80 
64.16 
26.55 
3.98 
-3.54 
3.98 
26.55 
64.16 
1 16.80 
94.25 
49.1 2 
19.03 
3.98 
3.98 
19.03 
49.12 
94.25 
86.73 
49.1 2 
26.55 
19.03 
26.55 
49.1 2 
86.73 
94.25 
64.16 
49.1 2 
49.12 
64.16 


312.70 

1251 .00 

1095.00 

938.10 

781.80 

625.40 
-629.10 
-786.90 
-944.80 

-1103.00 

-314.60 

-472.30 

-630.10 

-788.00 

-946.00 

0 

-157.40 

-315.10 

-472.90 

-630.80 

-788.80 

314.60 

157.40 
0 

-157.70 
-315.50 
-473.40 
-631 .40 

629.10 

472.30 

315.10 

157.70 
0 

-157.80 

-315.70 

-473.70 

786.90 

630.10 

472.90 

315.50 

157.80 
0 

-157.90 

944.80 

788.00 

630.80 

473.40 

315.70 

157.90 

1103.00 

946.00 

788.80 

631.40 

473.70 


-1083.00 

0 

-270.80 
-541 .60 
-812.50 
-1083.00 

910.20 

636.70 

363.10 

89.70 

910.70 

637.30 

363.80 
90.37 

-182.90 

910.80 

637.50 

364.20 
90.81 

-182.40 

-455.40 

910.70 

637.50 

364.30 
91.03 

-182.1 0 
-455.10 
-727.80 

910.20 

637.30 

364.20 
91 .03 

-182.10 

-455.00 

-727.60 

-999.90 

636.70 

363.80 
90.81 

-182.10 
-455.00 
-727.60 
-999.70 

363.10 
90.37 
-182.40 
-455.1 0 
-727.60 
-999.80 

89.70 
-182.90 
-455.40 
-727.80 
-999.90 













TABLE 11.- THERMAL LOADING ON RADIOMETER ANTENNA 
[Values are given relative to 80°F reference temperature] 


Temperature, °F 


5.1 

8.9 

8.4 

9.8 
21 .8 
17.2 
11 .9 

13.1 

14.4 

17.1 
16.6 

17.1 

13.4 

12.9 

15.9 

15.6 

18.5 

18.2 

17.2 

13.4 

14.5 

14.4 
14.1 

18.5 

19.7 

19.6 
31 .3 
13.5 

12.9 

12.7 

17.3 

17.0 

18.0 

19.1 

25.3 

6.9 

13.1 


Temperature, °F 


Temperature, °F 








TABLE 12.- SUMMARY OF ANTENNA REFLECTOR 
SENSITIVITY ANALYSIS 

109 nodes; 420 elements; 1 load condition; 
3 design variables 


(a) Solution time comparison 


Method 

Solution time, sec, for 
derivatives of - 

Displacement 

Stress 

Finite difference 

155 

160 

Semianalytical 

95 

114 


(b) Comparison of derivatives of 
center deflection 


Independent variable. 

Value of derivative of 

cross-sectional area 

center deflection 

Upper surface 

2.4 x 10 -4 

Diagonals 

-8.3 x 10- 5 

Lower surface 

-1 .8 x 1 0 -4 


63 





<T> 


EAL processors 



(one in 

central 

memory 

at a 

time) 

TAB 

K 

EQNF 

RMK 

AUS 

U4/VU 

ELD 

M 

SSOL 

DR 

DCU 

FSM 

E 

KG 

GSF 

SYN 

VPRT 

PS 

EKS 

LSK 

ES 

STRP 

PLTA 

Ell 

SEQ 

LSU 

PSF 

PXY 

PLTB 

PRTE 

TOPO 

RSI 

EIG 

PR 

SSBT 

DRSI 

TAN 

INV 

E4 

U1 

U3/RP2 


i 



Working storage area 


Data-handl ing utilities 



Data base complex 
(on auxiliary storage) 



Library 1 
Data set 1 


Library 2 
Data set 1 


Library 20 
Data set 1 



N 


N 


N 



Runstream 1 


Runstream 1 


Runstream 1 



M 


M 


M 


FI LE=L001 

FILE=L002 


FILE=L020 


*XQT TAB 
UPDATE= 1 
(TAB DATA) 

*CALL (l.CHNG DV) 

*XQT ELD 

(ELD DATA) 

*XQT K 
*XQT INV 


(b) Sample input file. 


*( l.CHNG DV) ENDCDV 
! I=NUMEL 
*LABEL 1 

*JLZ,-1(NUMEL,2)$L00P OVER ELEMENT NO.'S 
! N I = I - NUMEL 

$SELECT NEW VALUE OF DESIGN VARIABLE 
!DV=DS,1,"NIM(1,DV,DFN, MASK, MASK) 
$UPDATE DESIGN VARIABLE TO NEW VALUE 
{MEMBRANE THICKNESSES 
SA 

"NI" "DV" 

♦JUMP 1 

♦LABEL 2 

♦RETURN 
♦ENDCDV 


(a) System configuration. 


(c) Sample runstream data set. 


Figure 1.- The EAL system terminology and configuration 





(a) Finite-element model of cantilevered 
channel-section beam* 



(b) Section A-A showing cross-sectional 
geometry and structural variables* 

Figure 2.- Cantilevered channel-section beam used to illustrate analytical 

derivative method. 




View along 
negative X-axis 
(not drawn to scale) 


Figure 3.- Geometry and node numbering for swept wing. All dimensions are in inches. 


CO 





Figure 4.- Element numbering for swept wing. 
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Y 



X 


(a) Overall model. 

Figure 8.- Finite-element model of radiometer-antenna reflector. 
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1 


6 


12 


19 


27 



(b) Upper surface. 
Figure 8.- Continued. 
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66 


71 


77 


84 



(c) Lower surface. 
Figure 8.- Concluded. 
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Start 



Figure 9.- Flowchart for system of runstreams to calculate structural-sensitivity 
derivatives with semi analytical method. Numbers indicate appropriate sections 
of runstreams in appendix. 







Figure 10.- Flowchart of section of runstream system which calculates 
derivatives of stiffness, mass, and geometric stiffness matrices. 
Numbers indicate appropriate sections of runstreams in appendix. 


75 




Figure 11.- Flowchart of section of runstream system which details 
calculation of displacement and stress derivatives. Numbers 
indicate appropriate sections of runstreams in appendix. 





Figure 12.- Flowchart of section of runstream for calculation of 
derivatives of vibration and of buckling eigenvectors and 
eigenvalues. Numbers indicate appropriate sections of 
runstreams in appendix. 


77 



1. Report No. 

NASA TM- 85743 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 

STATIC AND DYNAMIC STRUCTURAL-SENSITIVITY DERIVATIVE 
CALCULATIONS IN THE FINITE-ELEMENT-BASED ENGINEERING 
ANALYSIS LANGUAGE (EAL) SYSTEM 


5. Report Date 
March 1984 


6. Performing Organization Code 

506-53-53-07 


7. Author(s) 

Charles J. Camarda and Howard M. Adelman 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


8. Performing Organization Report No. 
L- 1 5659 


10. Work Unit No. 


11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 


This paper describes the implementation of static and dynamic structural-sensitivity 
derivative calculations in a general purpose, finite-element computer program denoted 
the Engineering Analysis Language (EAL) System. Derivatives are calculated with 
respect to structural parameters, specifically, member sectional properties including 
thicknesses, cross-sectional areas, and moments of inertia. Derivatives are obtained 
for displacements, stresses, vibration frequencies and mode shapes, and buckling 
loads and mode shapes. Three methods for calculating derivatives are implemented 
(analytical, semianalytical, and finite differences), and comparisons of computer 
time and accuracy are made. Results are presented for four examples: a swept wing, 


a box beam, a stiffened cylinder with a cutout, and a space radiometer-antenna truss. 


17. Key Words (Suggested by Author(s)) 


18. Distribution Statement 


Structural sensitivity analysis 
Large space antenna 
Structural optimization 
Finite-element method 
Thermal loads 

Engineering Analysis Language 

Unclassified - Unlimited 

Subject Category 39 

19. Security Classif. (of this report) 

I 20. Security Classif. (of this page) 

21. No. of Pages 

22. Price 

Unclassified 

Unclassified 


78 

AOS 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley , 1984 





















National Aeronautics and 
Space Administration 


THIRD-CLASS BULK RATE 


Washington, D.C. 
20546 


Official Business 

Penalty for Private Use, $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 





POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 




